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Z.  IMTRODUCTIOM 


The  objective  of  this  thesis  xs  to  design,  construct  and 
calibrate  a  pair  of  reduced  order  force-torque  transducers 
based  on  force  sensing  resistor  technology.  This  is  one 
branch  of  a  two  branch  robotics  research  project.  The  work 
done  in  parallel  is  a  separate  study  of  how  a  hybrid  force 
override  rate  control  system  can  be  used  to  enable  operation 
of  robotic  manipulators. 

The  motivation  for  this  research  is  that  robotic 
applications  such  as  the  Space  Shuttle  Manipulator  Arm  require 
more  than  just  rate  and  position  control  to  perform  certain 
assembly  tasks.  The  ability  to  detect  and  control  interaction 
forces  and  moments  is  necessary  in  poorly  defined  environments 
such  as  space.  To  accomplish  various  tasks  without  damage  to 
the  surrounding  objects  or  the  manipulator,  the  contact  forces 
must  be  sensed  and  used  in  a  feedback  loop  to  regulate  the 
conventional  rate  and  position  control.  Further,  with  the  use 
of  force  override  control  systems,  human  intervention,  which 
may  be  hazardous,  impractical  or  even  impossible,  could  be 
entirely  eliminated.  This  offers  the  potential  for  reducing 
costs,  protecting  human  life  and  enabling  the  accomplishment 
of  tasks  never  before  possible. 
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In  order  to  develop  a  practical,  simple  and  durable 
system,  the  design  of  the  force-torque  transducer  in  this 
project  is  based  on  the  Force  flensing  Resistor^  (FSR)  rather 
than  the  strain  gauge.  FSRs  provide  a  relatively  ine>cpensive, 
vibration  insensitive,  durable  and  adaptable  source  of 
resistance  changes  from  various  force  applications.  To  harness 
this  technology  in  the  design  of  a  force-torque  transducer,  it 
was  necessary  to  determine  a  basic  configuration  which  enables 
the  desired  forces  and  moments  to  be  extracted  from  various 
sensor  signals.  A  design  which  enables  the  resolution  of  three 
forces  and  three  moments  was  developed.  This  design  was  then 
modified  to  provide  three  forces  and  only  one  moment  in  the 
likely  scenario  of  applications  requiring  only  these  four 
commonly  used  degrees  of  freedom.  Also,  the  redundancy  of  the 
basic  six  degrees  of  freedom  transducer  was  analyzed  keeping 
in  mind  the  practical  use  of  such  a  device.  A  scheme  allowing 
multiple  levels  of  redundancy,  in  the  event  of  individual 
sensor  failure,  was  developed  to  demonstrate  the  functional 
reliability  of  a  force-torque  transducer  designed  using  force 
sensing  resistor  technology. 

Finally,  the  four  degree  of  freedom  transducer  and 
associated  electronic  interfaces  were  constructed  and  tested 
in  two  versions.  One  was  a  design  requiring  no  alterations  to 
the  readily  available  FSR.  The  other  required  a  new,  shorter 
FSR  not  currently  available  from  the  manufacturer.  To  overcome 
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this  obstacle  the  available  FSRs  were  modified  in  order  to 
build  a  force-torque  transducer  of  reduced  proportions.  This 
was  done  due  to  the  payload  limitations  of  the  PUMA  560 
robotic  manipulator  arm. 


ZZ.  PRELZMZIOAy  WORK 


A.  FORCE  OVBRRZDB  RATE  COHTROL 

1 .  Coiio«pt 

The  general  concept  of  a  force  override  rate  control 
system  is  that  the  contact  forces  between  a  robotic 
manipulator  end  effector  and  the  object  being  moved  can  be 
used  to  generate  a  feedback  signal.  This  signal  is  in  the 
form  of  a  voltage  change  which  is  used  to  alter  the  control 
input  signal  of  a  conventional  rate  or  position  control 
system.  This  gives  conventionally  controlled  manipulators  the 
ability  to  reflect  the  forces  applied  by  the  operator  on  the 
control  device. 

2 .  Background 

This  concept  has  been  tested  using  a  single  degree  of 
freedom  hydraulic  system  [Ref.l:pp.  5-12],  a  single  degree  of 
freedom  PUMA  560  manipulator  system  and  a  three  degree  of 
freedom  PUMA  560  manipulator  system  [Ref.  2;pp.  27-54] .  Figure 
1  provides  a  general  description  of  the  hydraulic  system  [Ref. 
1 : p . 6 ] .  Figure  2  shows  the  general  arrangement  of  the  PUMA 
system  [Ref.  3:p.  23].  These  systems  utilize  resistance  strain 
gauges  to  generate  voltage  signals  which  are  proportional  to 
the  end  effector  contact  forces. 
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Figur*  2 


on*  DOF  For  cm  Ovmrridm  R«t*  Control  of  *  PUMA  560 
Manipulator  [R*f.  2:p.  34] 
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In  general  a  strain  gauge  is  used  to  deternine  the 
actual  stress  experienced  on  the  surface  of  the  member  to 
which  it  is  bonded.  The  device  consists  of  a  folded  wire 
assembly  that  exhibits  a  resistance  change  as  the  length  and 
cross  sectional  area  change  according  to: 

R*pL/A  (1) 

Where  p  is  the  resistivity  of  the  material.  L  and  A  represent 
the  wire  length  and  cross  sectional  area.  Figure  3  shows  an 
enlarged  view  of  a  strain  gage.  These  resistance  changes  are 
small  and  temperature  effects  must  be  compensated  for  by  use 
of  a  second  unstrained  gauge  as  part  of  the  system.  Also,  the 


gauge  must  be  bonded  to  the  stressed  member  in  such  a  way  that 
all  the  strain  is  transferred  to  and  detected  by  the  gauge. 
Strain  gauge  rosettes  consisting  of  two  or  more  strain  gauges 
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with  different  orientations  are  used  when  stresses  are 
expected  in  more  than  one  direction.  Hooke's  law  is  used  to 
calculate  the  stresses  from  the  strains  which  are  proportional 
to  the  applied  loads. 

Force-torque  transducers  based  on  strain  gauge 
technology  have  been  applied  in  the  design  of  wrist  force 
transducers.  However,  these  tend  to  be  expensive  and  complex 
requiring  the  part  on  which  the  strain  gauge  is  mounted  to  be 
capable  of  deflection  [Ref.  l:p.  15].  Also,  such  systems  have 
a  high  susceptibility  to  noise  since  the  small  resistance 
changes  require  very  high  amplifier  gains  in  order  to  produce 
a  useful  voltage  signal.  Hence,  an  alternative  source  of 
voltage  output  proportional  to  force  input  was  sought  which 
was  low  cost,  robust  and  simple  to  construct,  calibrate  and 
use. 

B.  SENSOR  DEVELOPMENT 

1.  Force  Sensing  Resistor 
a.  Description 

(1)  Construction 

A  Force  Sensing  Resistor™  is  a  device 
manufactured  by  Interlink  Electronics  of  Santa  Barbara,  CA. 
Figure  4  shows  the  construction  of  a  typical  FSR.  A  sheet  of 
polymer  with  a  layer  of  semiconduct ive  ink  is  faced  with 
another  sheet  of  polymer  overlaid  with  a  conducting  pattern 
consisting  of  a  set  of  interconnecting  electrodes.  These 
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Figur*  4 

Typical  Force  Sanslng  Resistor  [Ref.  4:p.  1] 

electrodes  are  shunted  by  the  semiconducting  ink  when  the 
layers  are  brought  together . [Ref .  4:pp.  1-2] 

(2)  Characteristics 

The  FSR  has  a  very  high  no  load  resistance  on 
the  order  of  1  Mohm  or  more.  As  the  layered  device  is 
compressed  the  resistance  drops  proportionally.  Figure  5  shows 
a  typical  plot  of  resistance  verses  force  which  is 
significantly  more  sensitive  than  the  output  of  a  strain 
gauge.  Thus,  the  susceptibility  to  noise  is  immensely 
reduced. [Ref .  4:pp.  1-2] 

Figure  6  displays  the  results  of  a  longevity 
test  conducted  by  the  manufacturer  demonstrating  the 
durability  of  a  typical  FSR.  Interlink  Electronics  also 
reports  that  the  standard  temperature  range  is  up  to  about 


170”C.  Sensors  capable  of  withstanding  higher  temperatures 
can  be  manufactured.  FSRs  are  also  relatively  insensitive  to 
humidity . [Ref .  4:p.  3] 


FSR  Forcc/Rtsisunc*  Characteristic 


Figure  S 

Typical  F8R  Force/Resistanee  Characteristic  [Raf.4tp.  1] 


FSR  Longawity  Teat 


Tm  ConoMdi:  MT  Aa.  oreulv  praoa: 

10  man  sBikM  ■  52  9B. 
03  MC/npatton 


Figure  6 

FSR  Longevity  Test  [Ref.4:p.  3] 
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FSRs  can  be  readily  manufactured  in  various 
sizes  as  single  sensors  or  arranged  in  arrays.  It  is  also 
possible  to  vary  the  force  range  and  the  resistance  range  to 
meet  specific  requirements  during  fabrication.  Interlink 
Electronics  suggests  that  it  is  best  to  use  the  mechanical 
design  to  establish  the  force  within  a  useable  range. [Ref. 
4:p.  2] 

An  FSR's  performance  is  somewhere  between  a 
true  force  sensor  and  a  pressure  sensor.  For  a  given  force, 
a  pressure  sensor  yields  an  output  inversely  proportional  to 
the  area  of  the  applied  force.  Whereas  a  force  sensor  will 
yield  a  constant  output  regardless  of  the  area  and 
distribution  over  which  the  given  force  is  applied.  If  the 
force  distribution  area  is  smaller  than  the  FSR  active  area 
but  large  compared  to  the  spacing  between  the  interconnecting 
electrodes,  then  the  resistance  output  will  vary  approximately 
as  the  reciprocal  of  the  square  root  of  the  applied  force 
area.  However,  when  the  area  of  the  applied  load  is  larger 
than  the  FSR  active  area  the  device  can  -  be  used  as  a 
pressure  sensor. [Ref.  4:p.  2] 

Sensor  electrical  interfaces  can  be  simple 
because  typical  resistance  changes  are  relatively  large; 
therefore,  these  circuits  do  not  require  a  resistance  bridge 
as  do  strain  gauge  circuits.  Also,  the  impedance  is  almost 
entirely  resistive.  Depending  on  the  circuit  used,  sensor 
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response  signals  are  typically  on  the  order  of  0-10 
volts. [Ref .4:p.  4] 

Force  sensing  resistor  technology  offers  many 
advantages  over  conventional  strain  gauge  technology.  In  an 
attempt  to  develop  the  force  override  rate  control  system  as 
a  low  cost,  rugged  and  credible  concept,  a  suitable  force- 
torque  transducer  must  be  developed.  Hence,  the 
implementation  of  readily  available  FSRs  is  a  logical  choice 
for  a  simple,  robust  and  inexpensive  transducer  design. 

2 .  Transducer 
a.  Principle 

The  use  of  force  sensing  resistor  technology  was 
suggested  to  the  project  sponsor  at  NASA.  Figure  7  shows  the 
sponsor's  suggested  conceptual  design  for  a  suitable  force- 
torque  transducer  using  readily  available  FSRs  [Ref.  l:p.  16]. 
Initially,  it  was  proposed  that  the  sensors  be  mounted  on  an 
inner  cube  which  would  be  surrounded  with  a  semi -compliant 
material  like  RTV  and  encased  in  a  larger  cube.  When  the 
joystick  is  moved,  some  of  the  FSRs  will  experience  local 
normal  forces  and  thus  a  change  in  resistance.  When  the 
appropriate  electronic  interface  is  used,  a  corresponding 
voltage  change  will  be  produced  resulting  in  sensor  signals 
which  reflect  the  applied  forces  and  moments.  Ultimately,  the 
sensors  would  be  located  in  such  a  way  that  the  applied  forces 
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and  moments  could  be  detected  and  decoupled  for  use  in  a  force 
override  rate  control  system. [Ref.  3:p.  11] 


Figure  7 

FSR-Based  Transducer  Conceptual  Design  [Ref.  l:p.  17] 
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b.  Sensor  Equations 

The  original  Biathenatically  generalized  description 
of  the  problem  was: 


In  this  equation  the  sensor  vector,  8,  represents  the  voltage 
output  for  each  sensor  when  the  force  and  moment  vector,  F,  ic 
applied,  and  n  represents  the  number  of  sensors  used  which 
depends  on  the  selected  configuration.  The  positive  and 
negative  directions  for  each  force  and  moment  component  had  to 
be  accounted  for  because  the  sensors  were  only  used  to  detect 
applied  compressive  loads.  Hence,  the  objective  was  to 
determine  the  right  number  and  location  of  sensors  and  the 
associated  A  matrix  in  order  to  resolve  the  12x1  force  and 
moment  vector  using: 

F=A'8  (3) 
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The  elements  of  the  A  matrix  would  have  to  be  determined  using 
a  calibration  procedure  which  measured  each  sensor  signal 
after  applying  known  forces  and  moments  individually.  [Ref. 
l:p.  18] 

Since  the  number  of  sensors  could  be  nore  than 
twelve,  the  A  matrix  may  not  be  invertible.  In  which  case  the 
pseudoinverse  must  be  utilized  according  to: 

F=(A^A) ‘A^S  (4) 

If  the  square  matrix  A’^A  is  of  full  rank  (i.e.  a  rank  equal  to 
the  length  of  the  force  and  moment  vector) ,  it  may  be  inverted 
and  used  as  shown  in  Equation  4  to  identify  the  forces  and 
moments.  [Ref.  l:p.  19] 

The  initial  model  of  a  control  device  to  detect  the 
forces  and  moments  was  that  of  a  planar  joystick  shown  in 
Figure  8.  An  equation  was  developed  for  each  sensor  using  the 
principle  of  superposition  and  modeling  each  FSR  as  a  linear 
spring.  The  principle  of  superposition  was  used  because  the 
application  of  a  force,  P,  at  the  top  of  the  joystick  would 
cause  the  device  to  translate  and  rotate.  Hence,  the 
following  equations  account  for  the  contribution  of  both 
effects  to  the  normal  force  on  each  spring: 


F,=P/4(H-d/a) 

(5) 

F2=P/4(l-d/a) 

(6) 

F3=-P/4(l+d/a) 

(7) 

F4=-P/4(l-d/a) 

(8) 
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Figur*  8 

Planar  Joyatiok  (Raf.  3:p.20] 


By  applying  a  pure  moment  about  the  y  axis  the  equations 
become : 

F,=F4=M/4a  (9) 

F2=F3=-M/4a  (10) 

The  sensor  coefficient  matrix  A  can  be  determined  for  a 
variety  of  three  dimensional  configurations  by  extending  this 
methodology  to  a  cube  and  determining  the  effects  of  each 
applied  force  and  moment.  [Ref.  l:pp. 17-21] 


c.  Sensor  Placement 

Since  the  A  matrix  is  dependent  upon  the  location 
of  the  sensors  on  the  cube,  the  first  task  in  this  iterative 
design  process  was  to  pick  a  number  and  distribution  of  FSRs 
on  the  cube.  Then  by  applying  the  forces  and  moments  of 
interest,  the  components  of  the  A  matrix  were  determined  using 
the  principle  of  superposition  outlined  earlier.  The  A^A 
matrix  was  then  computed  and  its  rank  checked.  If  the  results 
indicate  a  rank  equal  to  the  length  of  the  force  and  moment 
vector,  then  that  sensor  configuration  was  considered  capable 
of  being  used  to  resolve  the  forces  and  moments.  If  the  rank 
was  greater  than  the  vector  length,  then  too  many  sensors  were 
being  considered.  If  the  rank  was  less  than  the  vector 
length,  then  too  few  sensors  were  being  considered.  (Ref.  3;p. 
14] 

d.  Early  Results 

The  iterative  process  for  determining  the  right 
number  and  location  of  sensors  was  carried  out  for  a  full 
order  transducer  consisting  of  6  forces  and  6  moments.  This 
accounted  for  three  dimensional  forces  and  moments  with 
sensors  only  capable  of  registering  a  compressive  load.  A 
possible  placement  pattern  which  resulted  in  an  A^A  matrix 
with  a  rank  of  12  is  shown  in  Figure  9.  [Ref.  l;pp.  29-30] 

The  iterative  process  was  then  carried  out  in  a 
similar  fashion  for  a  reduced  order  transducer  capable  of 
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resolving  the  three  dinensional  forces  as  before,  but  a  moment 
only  about  the  z  axis.  This  was  done  because  not  all  robotic 


Figure  9 

Arrangement  for  Pull  order  Transducer  30] 

manipulator  tasks  require  moments  in  all  three  directions.  It 
is  common  to  only  require  a  torque  about  a  single  axis  in 
addition  to  the  forces  along  all  three  directions.  Again, 
since  the  FSRs  are  only  responsive  to  compressive  loads,  a 
separate  sensor  had  to  be  allowed  for  the  negative  and 
positive  directions  of  each  of  these  components.  The 
potential  sensor  layout  with  a  rank  of  8  is  shown  in  Figure 
10.  [Ref.  l:pp.  31-33] 
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It  Is  apparent  that  in  order  to  reduce  the  nunber 
of  sensors  required,  a  design  modification  is  necessary  which 
addresses  the  response  of  FSRs  only  to  compressive  loads. 


Figure  10 

Arrangement  for  Reduced  Order  Transducer  [Ref.  l:p.33] 

This  could,  in  effect,  halve  the  force  and  moment  vector 
length  by  not  requiring  sensors  for  negative  direction  forces 
and  moments. 
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III.  PRE-LOADED  8E1I80R  TRAM8DUCER  DEVELOPMEMT 


A.  8ZHQLE  F8R  PROTOTYPE 

1.  Obj«otiv« 

In  an  effort  to  resolve  the  problem  of  the  FSR's 
inability  to  register  a  signal  when  a  tensile  load  is  applied, 
the  concept  of  pre-loading  the  sensor  was  employed.  This 
required  the  sensor  to  be  loaded  compressively  with  some 
initial  force.  The  change  in  resistance  caused  by  adding  or 
removing  normal  force  applied  to  the  FSR  is  used  to  generate 
a  proportional  change  in  voltage.  Provided  that  the  initial 
pre-load  is  sufficient,  a  useable  range  is  created  in  which  a 
response  is  generated  to  both  tensile  and  compressive  force 
application. 

2.  Design  end  Construction 

In  order  to  test  this  concept  a  device  had  to  be 
created  which  would  show  a  reasonable  range  of  resistance 
change  for  varying  compressive  and  tensile-  normal  force 
application. 

a.  Mechanical  Device 

Figure  11  shows  the  prototype  used  to  test  the  pre¬ 
load  concept.  The  device  was  constructed  using  standard  .25” 
aluminum  plate.  As  shown  .25”  of  neoprene  material  was  used 
to  transmit  the  applied  load  to  the  FSR  which  was  mounted 
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using  Scotch™  467MP  High  Performance  adhesive.  The  FSR  was 
.25”  square,  product  #30-301,  readily  available  from  Interlink 
Electronics™.  A  standard  1.5”  10/32  hex-head  machine  bolt 
with  a  locking  nut  was  to  apply  the  initial  pre-load.  A 
detailed  drawing  is  provided  in  Appendix  A. 


Figure  11 

Single  Pre-loeded  FSR  Prototype 
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b,  Blsctronic  InterfacB 


In  order  to  obtain  a  voltage  change  proportional  to 
the  applied  load,  the  electronic  interface  pictured  in  Figure 
12  was  used.  The  circuit  consisted  of  a  -flO  volt  reference 
source,  various  values  of  and  an  LM324  op-  ^p.  The  output 
voltage  can  be  %nritten  as: 

Voitt“Vh£p/2«  t  l+Ro/Rps*]  (11) 

The  expected  voltage  swing  is  from  V|t£p.  [Ref.  5:p.  8] 


Figure  12 

F8R  Current-to-Voltage  Converter  [Ref.  5:p.  8] 
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3  .  R«spoilS«  T«St 

a.  Data  Collection 

The  prototype  was  configured  to  sit  on  a  table  top 
while  a  compressive  load  was  applied  by  adding  weights  to  a 
tray  on  top  of  the  pre-load  set  bolt.  Tensile  loads  were 
applied  by  turning  the  device  over  and  suspending  it  by 
supporting  the  base  on  two  cross  members  across  a  sufficient 
gap  between  two  adjacent  and  parallel  tables.  The  weights 
were  then  added  to  a  tray  suspended  from  the  pre-load  set 
bolt. 

Initially,  the  masses,  ranging  from  4.17  to  -1.97 
Ibm,  were  measured  in  English  units  and  the  corresponding 
voltages  were  read  using  the  OMEGA  BENCH™  data  acquisition 
system  and  APPLE  MACINTOSH™  microcomputer  configured  as  a 
voltmeter.  The  first  test  used  4.7  kohms  for  the  value  of  Ro 
in  an  effort  to  produce  a  curve  similar  to  that  shown  in 
Figure  13.  This  was  the  response  for  a  standard  FSB  using  a 
similar  interface  obtained  by  Interlink  Electronics™.  [Ref. 
5:p.  7] 

In  order  to  more  fully  understand  the  nature  of 
pre-loaded  FSRs,  similar  tests  were  conducted  using  the 
following  values  for  R^  in  kohms:  2.5,5.0,7.5,10.0,12.5.  The 
masses  used  during  these  tests  were  measured  in  metric  units 
for  convenience. 
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rlgure  13 

Force  vs.  Vot^r  Curves  [Ref.  5:p.  7] 

b.  Results 

The  data  taken  and  the  corresponding  plots  of 
applied  force  versus  voltage  response  are  provided  in  Appendix 
A.  The  plot  resulting  from  the  first  test  is  sho%m  in  Figure 
14.  The  curve  resembles  the  expected  results  shown  in  Figure 
13.  The  plot  clearly  shows  two  nearly  linear  regions 
separated  by  the  initial  no  load  set  points.  In  subsequent 
tests  the  initial  set  point  was  physically  unchanged. 
However,  different  values  for  Rq  resulted  in  curve  shifts. 
This  is  seen  in  the  variation  of  the  no  load  voltage  levels  in 
accordance  with  the  various  values  for  Rg.  Also,  depending  on 
the  value  of  the  initial  voltage,  the  useable  range  of  the 
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pre-loaded  sensor  was  affected.  Too  high  a  value  for  Rq 
caused  positive  saturation  as  in  the  12.5  kohm  case.  The 
conclusion  drawn  from  these  tests  was  that  the  FSRs  configured 
in  a  simple  pre-loaded  manner  with  a  simple  electronic 
interface  are  capable  of  registering  a  potentially  useable 
voltage  change  which  over  limited  regions  varies  almost 
linearly  with  applied  loads  in  both  the  compressive  and 
tensile  directions. 

B.  INITIAL  TRANSDUCER  CONTIODRATION 
1 .  Arrangement 

Since  each  FSR  can  now  be  used  to  detect  normal  force 
loading  in  both  the  compressive  direction  and  the  tensile 
direction,  the  12x1  force  and  moment  vector  in  Equation  2  is 
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reduced  to  a  6x1  vector.  The  mathematically  generalized 


description  of  the  problem  is  now: 
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Where  8  represents  the  sensor  voltage  output  and  F  is  the 
force  and  moment  column  vector.  The  guest  is  still  to 
determine  a  number  and  configuration  of  sensors  such  that  the 
6x1  force  and  moment  vector  can  be  resolved  by  inverting  the 
A  matrix  according  to: 

FssA**8  (13) 

If  the  sensor  configuration  is  such  that  only  6  sensors  are 
used,  then  the  A  matrix  will  be  square  and  the  pseudoinverse 
will  not  be  required  as  shown  previously  in  Equation  4. 

Figure  15  shows  a  three  dimensional  transducer  built 
upon  the  planar  joystick  concept  discussed  in  the  previous 
chapter.  The  plan  was  to  develop  the  equations  necessary  to 
determine  the  A  matrix  with  all  the  sensors  shown.  Then, 
systematically  cut  out  the  signals  from  various  groups  of 
sensors  and  check  the  rank  of  the  corresponding  A  matrix.  If 
the  rank  was  greater  than  or  equal  to  six,  then  the 
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configuration  was  considered  capable  of  resolving  the  6 
components  of  the  force  and  moment  vector.  The  goal  was  to 
determine  a  configuration  which  uses  the  minimum  number  of 


sensors  and  has  a  corresponding  invertible  A  matrix. 


Figure  15 

Initial  Transducer  Configuration 


2 .  Sensor  Equations 

Using  the  principle  of  superposition  mentioned  earlier 
for  the  planar  joystick,  equations  were  developed  for  forces 
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applied  in  the  x,  y,  and  z  directions  as  well  as  the  moments 
about  each  axis.  These  forces  and  moments  were  considered  to 
have  been  applied  at  the  top  of  the  joystick  handle.  First, 
a  force  applied  in  the  positive  x-direction,  F,,  results  in  a 
translation  of  the  inner  cube  as  well  as  an  incidental  moment 
about  the  y  axis.  Since  there  are  eight  sensors  on  the  x 
faces,  the  translation  force  contribution  is  simply: 

F'=F,/8  (14) 

The  sign  of  this  force  equation  for  each  sensor  depends  on 
whether  the  sensor  is  on  the  positive  or  the  negative  x  face. 
The  incidental  moment  created  by  F,  causes  activation  of  16 
sensors,  eight  on  the  x  face^  and  eight  on  the  z  faces.  Hence, 
the  resulting  applied  load  due  to  this  incidental  moment  can 
be  determined  for  each  affected  sensor  according  to: 

16F"a=M  (15) 


where , 


M=F.d 


(16) 


hence , 


F"=F,d/16a  (17) 

Again,  the  sign  of  this  force  equation  for  each  sensor  depends 
on  its  location.  Since  this  force  is  a  result  of  the 
incidental  moment,  the  location  of  the  sensor  with  r.ispect  to 
the  axis  of  rotation,  the  y  axis  in  this  case,  determines  the 
sign  of  F”. 
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Finally,  with  the  signs  determined  correctly  for  each 
activated  sensor's  force  equations,  the  normal  force 
experienced  by  each  FSR  resulting  from  the  application  of  F, 
is  the  sum  of  the  translational  force  and  the  vorce  from  the 
incidental  rotation: 

F=F'+F''  (18) 

This  rationale  can  be  repeated  for  a  force  applied  in 
the  y  direction,  F^.  However,  a  force  applied  in  the  z 
direction,  F^,  would  not  result  in  an  a  rotation  of  the  inner 
cube.  Therefore,  the  resulting  normal  forces  on  the  affected 
sensors  would  only  be  from  the  translation  of  the  inner  cube 
with  no  incidental  moment. 

For  a  pure  torque  applied  at  the  end  of  the  joystick 
in  the  x  direction,  M,,  there  would  be  sixteen  sensors 
affected,  8  on  the  y  face  and  8  on  the  z  face.  The  resulting 
force  Taagnitude  for  each  activated  sensor  is: 

F=M,/16a  (19) 

The  sign  for  each  of  these  activated  sensor  equations  depends 
on  its  location  with  respect  to  the  axxs  of  rotation.  Pure 
moments  applied  about  the  y  axis  and  about  the  z  axis.  My  and 
Mj  respectively,  result  in  similar  force  equations  for  each 
affected  sensor. 

The  equations  for  each  sensor  activated  as  a  result  of 
the  application  of  the  individual  components  of  the  force  and 
moment  column  vector  make  up  the  A  matrix  in  Equation  12.  The 
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complete  equation  for  this  configuration,  showing  the 


individual  elements  of  the  A  matrix,  is  as  follows: 
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3 .  Sensor  Placement 

This  A  matrix,  resulting  from  the  configuration  shown 
in  Figure  15,  was  entered  into  the  software  program  MATLAB™ 


with  the  assumed  values  for  d  and  a  as  4.5  and  1.5, 
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respectively  with  units  of  length.  The  the  rank  of  the  full 
matrix  was  then  determined  to  be  six.  This  full  rank  meant 
that  if  all  sensors  were  used,  the  pseudoinverse  could  be  used 
to  resolve  the  components  of  the  force  and  moment  vector,  F. 

Next,  the  sensors  on  the  negative  x  face  were 
eliminated  by  setting  the  X  matrix  elements  in  corresponding 
rows  equal  to  zero.  In  this  case,  it  was  simulated  that 
sensors  labelled  -1,  -2,  -3,  and  -4  had  no  output.  Thus,  the 
elements  in  rows  13  through  16  of  matrix  X  were  set  equal  to 
zero.  The  resulting  matrix  had  a  rank  of  six.  Hence,  this 
configuration  could  be  used  to  resolve  the  components  of  F. 

This  process  was  repeated  in  a  similar  fashion  by 
first  eliminating  additional  sensors  from  the  y  face,  and  then 
from  the  z  face.  The  result  was  a  configuration  with  sensors 
only  on  the  positive  x,  y,  and  z  faces  yielding  an  X  matrix 
with  a  rank  of  six.  Since  this  would  involve  a  total  of  12 
sensors,  the  X  matrix  would  not  be  square.  However,  the 
pseudoinverse  could  be  used  to  resolve  the  forces  and  moments. 

In  an  effort  to  reduce  the  number  of  sensors  required 
even  further,  this  iterative  process  of  removing  sensors  from 
the  configuration  was  continued.  The  next  removal  was  that  of 
the  positive  z  face  sensors.  This  resulted  in  an  X  matrix 
rank  of  five  which  meant  this  sensor  configuration  could  not 
be  used.  Hence,  two  sensors  were  added  back  to  the  positive 
z  face.  These  were  the  diagonally  arranged  sensors  numbered 
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10  and  12  in  Figure  15.  The  resulting  A  natrix  yielded  a  rank 
of  six.  This  indicated  a  useable  conbination. 

Using  this  rationale,  two  sensors  were  also  removed  in 
diagonal  patterns  from  both  of  the  positive  x  and  y  faces. 
The  sensors  numbered  2,  3,  5  and  8  were  eliminated.  The 
resulting  rank  was  five  indicating  an  unusable  combination. 
Hence,  the  iterative  process  was  continued  to  find  a  useable 
combination  of  only  six  sensors. 

4 .  Results 

Experimenting  with  diagonal  combinations  lead  to  the 
useable  configuration  sho%m  in  Figure  16  after  only  three  more 
iterations.  This  is  not  the  only  possible  combination,  but  it 


Figure  16 

Initial  Minimum  Sensor  Transducer  Configuration 
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is  one  that  yields  a  square  A  matrix  which  is  directly 
invertible  to  determine  the  force  and  moment  vector,  F.  For 
fewer  than  six  sensors,  the  rules  of  matrix  multiplication 
prohibit  the  formation  of  the  pseudo inverse.  Hence,  no  fewer 
than  six  sensors  were  considered  in  any  configuration. 

C.  COMPARISON  WITH  PRELIMINARY  WORK 

1 .  Arrangement 

For  comparison  purposes,  an  unsuccessful  configuration 
analyzed  in  the  preliminary  work  for  a  full  order  transducer 
without  using  the  pre~ loaded  sensor  concept  was  reexamined 
using  pre-loaded  sensors.  This  arrangement  is  shown  in  Figure 
17  [Ref.  l:p.  25].  The  goal  was  to  determine  what 
modifications,  if  any,  had  to  be  made  to  create  a 
configuration  using  only  six  sensors  which  would  produce  an  A 
matrix  with  a  full  rank. 

2 .  Sensor  Equations 

Using  this  configuration  and  the  fact  that  the 
individual  FSRs  are  capable  of  measuring  tensile  and 
compressive  loads,  a  new  set  of  governing  force  equations  were 
developed.  These  equations  were  determined  using  the  principle 
of  superposition  analogous  to  the  methods  used  in  the  previous 
section  by  individual  application  of  each  component  of  the 
force  and  moment  vector,  P. 

Since  only  12  sensors  were  considered  in  this  initial 
configuration,  in  contrast  to  the  24  of  the  previous  section. 
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rigur*  17 

12  8Msor  Pr«liaiiiary  Work  Coafigurotion  [Rof.  X:p.  25] 


the  complete  equation  is  as  follows: 
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since  there  are  half  as  many  sensors  over  which  the  forces  are 
distributed,  the  resulting  elements  of  the  A  matrix  are  twice 
as  large  as  in  Equation  20.  This  has  absolutely  no  effect 
since  multiplying  the  A  matrix  by  a  scalar  will  not  change  its 
rank. 

3 .  Sensor  Plaoement 

Using  MATLAB''^  as  before,  the  full  A  matrix  was 
entered  and  its  rank  determined  to  be  six.  The  length 
dimensions,  d  and  a,  where  assumed  to  be  5.0  and  1.0, 
respectfully.  This  full  rank  indicated  that  the  configuration 
could  be  used  to  resolve  the  six  components  of  the  force  and 
moment  vector,  F.  In  an  attempt  to  optimize  the  number  and 
placement  of  sensors,  the  sensors  on  the  negative  x,  y,  and  z 
faces  were  eliminated,  leaving  active  sensors  numbered  1,  2, 
5,  6,  11  and  12.  The  resulting  configuration  yielded  an  A 
matrix  with  a  rank  of  only  five.  Since  this  is  a  rank 
deficient  condition,  another  attempt  was  made  by  eliminating 
the  sensors  on  only  the  negative  x  and  y  faces,  and  the 
sensors  on  the  positive  z  face.  This  configuration  left 
sensors  numbered  1,  2,  5,  6,  9  and  10  active.  The  rank  of  the 
resulting  matrix  was  five,  again  a  deficient  value,  indicating 
an  unusable  combination. 

4 .  Results 

It  was  intuitively  determined  that  in  order  to  develop 
a  combination  utilizing  the  sensors  on  only  the  positive  faces 
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or  only  the  negative  faces,  this  configuration  nay  need  a 
slight  alteration.  Hence,  it  was  abandoned  in  search  of  a 
similar  yet  useable  configuration. 

D.  MODZFIBO  PRBLIMZMARy  COMPiaDRATZOlf 
1 .  JUrraagement 

Considering  the  previous  configuration.  Figure  18 
displays  a  similar  arrangement.  The  only  difference  is  the 


r' 

z 

Figure  18 

Modified  12  Sensor  Prelimiaery  Work  Configuration 

orientation  of  the  z  face  sensors.  Diagonally  located  sensors 
numbered  11  and  12,  located  on  the  positive  z  face,  have  been 
rotated  90°  clockwise  if  viewed  looking  in  the  positive  z 
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direction.  Similarly,  sensors  numbered  9  and  10  have  been 
shifted  90°  clockwise  if  viewed  in  the  positive  z  direction. 

2 .  seasor  Equations 

In  this  case  the  sensor  equations  which  comprise  the 
elements  of  the  A  matrix  were  very  similar  to  those  of  the 
previous  section.  Since  the  number  of  sensors  was  the  same, 
the  only  change  was  in  the  signs  of  the  various  force 
contributions  as  a  result  of  the  changed  sensor  orientation. 
The  complete  equation  describing  the  new  situation  is: 
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As  expected  the  only  changes  were  the  signs  of  the  elements 
corresponding  to  the  application  of  F,  and  on  sensors 


numbered  9,  10,  11,  and  12. 
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3 .  8«iisor  PlacMiMt 

MATLAB^  was  used  to  first  check  the  rank  of  the  A 
matrix  with  all  12  sensors  activated.  This  value  was  six, 
again  indicating  a  useable  configuration.  Next,  the  sensors 
located  on  the  negative  x,  y,  and  z  faces  were  again 
eliminated.  The  rank  of  the  associated  matrix  was  determined 
to  be  six.  With  this  rank,  the  configuration  is  useable. 

4 .  Results 

A  rank  of  six  indicates  a  full  rank  allowing  the  A 
matrix  to  be  directly  Inverted  to  determine  the  forces  and 
moments  according  to  Equation  3.  Since  the  minimum  number  of 
sensors  possible  is  six,  this  configuration  was  determined  to 
be  one  of  many  acceptable  choices.  Using  the  trial  and  error 
methodology  it  is  possible  to  determine  other  configurations 
using  only  six  sensors  which  yield  an  A  matrix  capable  of 
direct  inversion. 
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IV.  RBDUXDMICT  MOOiYSIS 


A.  GOAL 

Given  a  transducer  with  four  sensors  mounted  on  each  face 
of  the  Inner  cube,  there  is  more  than  one  combination  of  the 
minimum  number  of  sensors,  six,  which  will  provide  an  A  matrix 
with  a  full  rank  used  to  resolve  the  six  components  of  the 
force  and  moment  column  vector,  V.  It  would  be  advantageous 
to  develop  a  system  that  would  automatically  find  another 
useable  combination  of  remaining  sensors  necessary  to  continue 
operation  in  the  event  of  a  sensor  failure  while  in  use. 
Hence,  the  aim  was  to  develop  an  algorithm  which  would  search 
for  and  identify  various  combinations  of  six  sensors  that 
could  be  used  to  identify  the  force  vector,  F.  This  was  to  be 
done  with  a  24  sensor  transducer  utilizing  some  initial 
useable  configuration  of  six  sensors  subject  to  random  FSR 
failures  which  would  drive  the  system  to  search  for  an 
alternate  configuration.  The  number  of  possible  combinations 
would  be  described  by  the  following  equation: 

C(n,r)*  nl/ (n-r) !rl  (23) 

Where  n  represents  the  total  number  of  sensors  and  r 
represents  the  number  of  sensors  in  each  group  necessary  to 
resolve  the  force  and  moment  vector.  [Ref.  6:p.  1-8] 
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In  this  case,  n  equals  24  and  r  equals  6.  After  doing 
the  calculation,  the  resulting  number  of  combinations  is 
134,596.  Potentially,  this  many  combinations  would  have  to  be 
analyzed  to  find  another  configuration  which  would  yield  a 
full  rank  A  matrix.  However,  as  sensors  fail  the  nximber  of 
possible  combinations  decreases.  It  was  decided  that  with  this 
finite  limit  to  the  number  of  possible  search  iterations,  a 
FORTRAN  code  would  be  developed  to  mimic  the  methodology  of 
the  previous  chapter  in  finding  a  useable  combination. 

B.  THEORY 

1.  Sensor  Configuration 

As  a  prelude  to  developing  this  FORTRAN  code,  a  24 
sensor  configuration  had  to  be  developed  with  a  convenient 
numbering  system.  Also,  it  was  decided  to  use  a  configuration 
that  could  be  initially  broken  into  four  separate,  six  sensor 
configurations  which  were  symmetrical  and  individually  capable 
of  resolving  the  three  components  of  force  and  three 
components  of  moment.  The  purpose  of  these  subsystems  was  to 
initially  provide  four  layers  of  redundancy.  In  the  event  of 
a  sensor  failure,  a  quick  switch  to  another  subsystem  would 
allow  continued  use  of  the  joystick  without  the  delay  of 
randomly  searching  for  other  possible  combinations  of  six 
working  sensors  from  the  remaining  sensors.  This  could  be 
repeated  until  no  intact  subsystems  remained. 
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To  do  this,  a  24  sensor  transducer,  similar  to  that 
shown  in  Figure  15,  was  used.  The  initial  numbering  system 
thought  to  be  best  was  a  simple  numbering  of  the  sensors  from 
#1  through  #24.  This  is  shown  on  the  inner  cube  faces  in 
Figure  19.  Since  it  was  shown  in  the  previous  chapter  that 


First  Sequential  Cube  Face  Humbering  Scheme 

only  half  the  sensors  on  the  positive  cube  faces  were  required 
to  develop  a  useable  configuration,  the  symmetrically 
analogous  configurations  were  obvious  choices  for  other 
possible  configurations.  Hence,  the  four  separate 
configurations  were  comprised  of  two  positive  face  systems  and 
two  negative  face  systems.  As  shown  in  Figure  19,  these 
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systens  were  respectively  labeled:  A-positive,  B-positive,  A- 
negative  and  B-negative. 

To  develop  a  numbering  system  such  that  the  complete 
24x6  A  matrix  could  be  conveniently  broken  into  the  four 
separate  6x6  smaller  matrices  each  with  a  rank  of  six  as  a 
result  of  the  four  separate  six  sensor  configurations,  a  trial 
and  error  procedure  was  necessary.  First,  the  numbers  of  the 
sensors  were  arranged  such  that  A*positive  system  had  sensors 
il-06;  the  A-negative  system  had  sensors  #7-*#12;  the  B- 
positive  system  had  sensors  #13-#18;  and  the  B>negative  system 
had  sensors  #19-#24  as  shown  in  Figure  19.  Then,  in  a  manner 
similar  to  that  in  the  previous  chapter,  the  equations 
representing  the  normal  forces  on  each  sensor  were  developed 
as  a  result  of  the  individual  application  to  the  components  of 
F.  As  in  the  previous  chapter,  these  equations  were  placed  in 
the  A  matrix  and  values  for  the  length  dimensions  d  and  a  were 
assumed  to  be  4.5  and  1.5,  respectively.  The  matrix  rank  was 
determined  using  MATLAB™  to  be  six.  Next,  the  four 
subsystems  were  formed  into  smaller  6x6  matrices  and  checked 
for  ranks  of  six.  However,  the  B-positive  and  the  B-negative 
subsystem  matrices  each  had  ranks  of  only  five.  Hence,  this 
configuration  was  not  capable  of  providing  the  four  levels  of 
redundancy  initially  sought. 

To  find  a  configuration  that  would  yield  four  useable 
subsystems,  the  sensor  niunbering  scheme  had  to  be  changed.  As 
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shown  in  Figure  20,  the  first  such  change  was  to  swap  #5  with 
#17,  #6  with  #18,  #11  with  #23  and  #12  with  #24.  This  schese 
failed  to  produce  the  four  matrices  with  ranks  of  six.  In  this 


Second  Sequential  cube  race  Nunbering  Scheme 

case  the  A-positive  and  the  A**negative  subsystems  were  rank 
deficient.  In  an  effort  to  find  other  arrangements  which 
consisted  of  four  useable  subsystems  sequentially  numbered, 
this  iterative  procedure  was  continued.  Table  1  summarizes 
the  results  of  all  the  sensor  number  switches. 

Thr  third  scheme,  shown  in  Figure  21,  was  similar  to 
the  first  scheme  except  that  the  switch  of  #6  with  #17  and  #12 
with  #23  was  made.  The  resulting  matrices  were  of  full  rank. 
Hence,  this  was  the  numbering  scheme  chosen  for  the  FORTRAN 
code  development. 
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TABLE  1 
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Successful  Sequential  Cube  Face  Numbering  scheme 
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2 .  Algoritla 


The  general  procedure  followed  to  provide  multiple 
layers  of  redundancy  for  a  force-torque  transducer  is  now 
summarized.  It  began  with  a  24  sensor  transducer  consisting 
of  a  2x2  array  of  FSRs  on  each  of  the  six  inner  cube  faces. 
The  procedures  of  the  previous  chapter  were  utilized  to 
determine  the  normal  forces  on  each  sensor  resulting  from  the 
individual  application  of  the  three  dimensional  forces  and 
moments.  Nominal  length  dimensions  were  assumed  for  d  and  a 
and  the  resulting  24x6  matrix  was  multiplied  by  a  scalar  value 
of  24  to  produce  convenient  elements  for  display  purposes. 

The  rows  of  this  large  matrix  were  sectioned  into  four 
invertible  6x6  matrices,  arranged  sequentially,  representing 
individual  subsystems  capable  of  being  used  to  resolve  the 
force  and  moment  column  vector,  F,  according  to  Equation  3. 
Each  row  represented  a  single  sensor,  and  these  subsystems 
were  arranged  sequentially  for  FORTRAN  programming 
convenience.  The  ranks  of  these  resulting  6x6  matrices  were 
determined  to  ensure  invertibility. 

A  sensor  was  then  randomly  selected  for  failure.  The 
algorithm  had  to  then  determine  which  subsystem  this  FSR  was 
in  and  switch  to  another  intact  subsystem.  This  step  was 
repeated  until  no  intact  subsystems  remained.  A  determination 
was  then  made  as  to  whether  or  not  there  remained  enough 
active  sensors  to  conduct  a  random  search  for  a  combination  of 
6  remaining  sensors  to  provide  a  6x6  invertible  matrix  for  use 


44 


in  resolving  F  as  discussed.  If  there  were  at  least  6 
remaining  sensors,  comprised  of  at  least  one  sensor  remrining 
in  each  the  x,  y,  and  z  directions,  the  random  search  was 
commenced.  This  search  consisted  of  randomly  picking  6  of  the 
remaining  sensors,  again  with  at  least  one  in  each  direction. 
The  matrix  rows  associated  with  these  FSRs  were  then 
sequentially  arranged  into  a  6x6  matrix.  The  rank  of  this 
matrix  was  calculated  in  order  to  determine  invertiblilty .  If 
the  rank  was  six,  then  the  algorithm  requested  to  fail  another 
FSR  and  conduct  another  search.  If  the  rank  was  Less  than 
six,  then  a  different  combination  of  six  of  the  remaining 
sensors  was  chosen  and  the  rank  was  checked.  This  was 
repeated  until  a  successful  6x6  matrix  was  produced  or  a 
sufficient  number  of  iterations  had  occurred  in  which  case  the 
algorithm  could  be  terminated. 

C.  FORTRAN  PROGRAM 

The  flow  chart  for  this  algorithm  is  presented  in  Figure 
22.  The  resulting  FORTRAN  code  was  written  and  is  presented 
in  Appendix  B.  It  was  based  on  the  numbering  arrangement 
displayed  in  Figure  21.  The  values  assumed  for  the  length 
dimensions  shown  Figure  15,  d  and  a,  were  4.5  and  1.5, 
respectively.  The  resulting  24x6  A  matrix  was  multiplied  by 
the  scalar  value  24  in  order  to  produce  a  simple  numerical 
format. 
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rigur*  22 

Redundancy  Analysis  Flowchart 


D.  RESULTS 


The  resulting  FORTRAN  program  ran  successfully.  In  most 
cases,  many  combinations  were  generated  in  a  matter  of 
seconds.  The  slowest  part  to  the  process  was  operator 
interaction  and  screen  display.  A  sample  output  is  provided 
after  the  code  listing  in  Appendix  B. 

The  output  shows  a  listing  of  the  original  24x6  A  matrix 
with  the  option  to  correct  or  change  any  elements.  Next  the 
sequentially  ordered  subsystems  and  their  corresponding  ranks 
are  shovm.  The  program  then  randomly  picks  and  eliminates  a 
sensor.  The  display  shows  the  sensor  chosen  and  its 
corresponding  subsystem.  The  output  then  shows  which 
subsystems  have  not  yet  been  disturbed.  This  is  repeated  at 
the  operators  discretion  until  no  more  undisturbed  subsystems 
remain.  At  that  point  the  display  shows  the  program  shifting 
to  the  mode  of  randomly  picking  six  of  the  remaining  FSRs  and 
generating  a  new  6x6  subsystem  matrix  with  a  rank  of  six.  The 
number  of  iterations  required  to  generate  this  subsystem, 
which  must  meet  specific  requirements,  is  also  shown.  The 
operator  is  given  the  opportunity  to  allow  the  program  to 
eliminate  another  FSR  and  repeat  this  random  generation  of  a 
useable  6x6  matrix  from  the  remaining  FSRs.  This  process 
continues  provided  there  are  at  least  six  sensors  remaining 
with  at  least  one  in  each  direction. 


47 


V 


RBDUCBO  OROBR  TRAM8DUCBR  THBCHtY 


A.  MOTZVATZOM 

Many  robotic  manipulator  tasks  require  only  single  axis 
moments.  An  example  would  be  the  act  of  inserting  and 
tightening  a  threaded  screw  or  bolt.  The  transducer 
configurations  considered  thus  far  have  been  for  three  forces 
and  three  moments  about  each  of  the  x,  y,  and  z  axes.  In 
order  to  simplify  the  transducer  design  and  construction,  an 
analogous  configuration  was  developed  for  a  reduced  order 
transducer  requiring  torques  about  only  the  z  axis  and  the 
three  dimensional  forces. 

B.  8BM80R  CORTZOURATZON 
1 .  Arrangement 

The  concept  of  a  reduced  order  transducer  was 
investigated  in  the  preliminary  work  which  considered  each  FSR 
capable  of  detecting  only  compressive  loads  [Ref.l:pp.  31-31]. 
However,  when  the  concept  of  the  pre-loaded  sensor  is 
utilized,  the  analysis  is  similar  to  that  used  for  the  full 
order  transducer  configuration  discussed  in  Chapter  III. 

Since  it  is  desired  to  detect  and  resolve  three  forces 
and  only  one  moment,  the  force  and  moment  column  vector,  Y, 
has  only  four  elements.  The  reduced  order  mathematical 
generalization  of  the  problem  is  therefore: 
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*1.1  *1.2  *1.2  *1.4 

*2.1  •  •  •  F, 

•  •  •  •  K  F, 

*m.l  *m.2  *%3  *b.4  (23) 

As  before,  8  represents  the  sensor  voltage  output  and  F  is  the 
force  and  moment  column  vector.  Also,  since  it  is  desired  to 
minimize  the  number  of  sensors  required  to  resolve  the 
components  of  the  4x1  P  vector,  n  will  be  four.  The  goal  was 
to  choose  a  configuration  which  will  yield  an  invertible 
square  A  matrix. 

Figure  18  showed  the  modified  preliminary  transducer 
configuration  which  was  successful  in  providing  an  A  matrix 
with  a  full  rank.  Hence,  in  order  to  obtain  a  reduced  order 
transducer,  the  iterative  process  of  removing  sensors  and 
their  corresponding  A  matrix  rows  began  with  this  layout. 

2 .  Equations 

As  described  in  Chapter  III,  the  principle  of 
superposition  was  employed  to  develop  the  equations  describing 
the  normal  forces  on  each  sensor  as  a  result  of  the 
application  of  the  three  dimensional  forces.  Also,  equations 
were  developed  to  reflect  the  application  of  moment  about  only 
the  z  axis.  The  resulting  complete  equation  is  similar  to 
Equation  22.  The  only  difference  is  that  the  two  A  matrix 
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columns  which  were  a  result  of  the  application  of  M,  and  My, 
are  now  gone.  The  new  equation  is  as  follows: 
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3 .  Sensor  Plaeement 

First,  the  values  for  d  and  a  were  assumed  to  be  5  and 
1,  respectively.  Then  the  complete  A  matrix  was  put  into 
MATLAB^*^  and  the  rank  was  calculated  to  be  four.  As  expected, 
this  indicated  a  sufficient  configuration  to  resolve  the 
components  of  F.  In  an  effort  to  reduce  the  number  of  sensors 
in  the  layout,  all  the  sensors  on  the  negative  faces  of  the 
cube  were  eliminated  by  setting  their  associated  rows  in  the 
A  matrix  to  zero.  This  resulted  in  a  configuration  with  six 
sensors  on  the  positive  face  of  the  cube  which  yielded  an  A 
matrix  with  a  rank  of  four. 
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The  layout  was  then  reduced  to  five  sensors  by  next 
removing  one  of  the  two  remaining  sensors  on  the  positive  z 
face.  The  choice  was  sensor  #11.  The  associated  K  matrix 
yielded  a  full  rank  of  four.  Next,  sensor  #12,  the  last 
remaining  z  face  sensor,  was  removed.  The  resulting  A  matrix 
was  rank  deficient  indicating  a  combination  incapable  of 
resolving  F.  Hence,  this  sensor  was  returned  to  the 
conf igurat ion . 

Sensor  #5,  ay  face  FSR,  was  the  next  choice.  The 
resulting  configuration  generated  an  A  matrix  with  a  rank  of 
four  capable  of  being  used  to  resolve  the  components  of  the 
force  and  moment  column  vector.  For  continuity,  this  sensor 
was  returned  and  sensor  #1,  a  positive  x  face  sensor,  was 
removed.  This  also  yielded  an  A  matrix  with  a  rank  of  four 
indicating  another  successful  combination.  Finally,  a 
configuration  which  was  nearly  a  mirror  image  to  the  layout 
missing  sensor  #5  was  used  generate  an  A  matrix  of  full  rank. 

4 .  Results 

After  the  eight  iterations  discussed  above,  two 
useable  sensor  configuration  were  determined.  These  layouts 
consist  of  the  minimum  number  of  pre-loaded  sensors,  four, 
necessary  to  resolve  the  elements  of  the  4x1  force  and  moment 
column  vector  using  Equation  13.  Figures  23  and  24  display 
the  resulting  transducer  configurations. 
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Figur*  23 

Poaitiv*  Fmoa  Pra*Loadad  Sanaor  Raduead  Ordar  Traasduoar 


Figura  24 

Nagativa  Faea  Pra-Loadad  Sansor  Raducad  Ordar  Traasduoar 
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VI.  PROTOTTPB  nUUISDUCBR 


B.  4  ZMCH  REDUCED  ORDER  VERSION 

1.  Obj«otlv« 

It  was  decided  to  build  a  suitable  prototype  force- 
torque  transducer  for  the  reduced  order  sensor  configuration. 
The  layout  to  be  used  is  pictured  in  Figure  23  and  was  sho%m 
theoretically  capable  of  being  used  to  resolve  forces  in  three 
directions  and  moments  about  a  single  axis.  After 
calibration,  a  physical  model  could  be  used  to  qualitatively 
demonstrate  the  capability  of  resolving  the  associated  force 
and  moment  column  vector  given  some  arbitrary  applied  loads  to 
the  joystick  handle. 

2.  Design  and  construction 
a.  MBchanical  Design 

Appendix  C  consists  of  the  detailed  machine 
drawings  of  a  4"  force-torque  transducer.  Figure  25  shows  a 
general  view  of  the  prototype  design.  The  outer  cube  was 
constructed  using  .25"  aluminum  plate  surrounding  a 
plexiglass,  smaller  solid  cube  with  an  attached  plexiglass 
joystick  handle.  The  handle  was  fastened  using  a  pin  and 
epoxy.  Phillips  head  4x1/2"  tapping  screws  were  used  to 
secure  the  walls  of  the  outer  cube. 
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Figur*  2S 

0«ii«ric  Prototyp*  Transdue«r 


The  4"  outside  dimension  was  chosen  in  order  to 
allow  for  a  3"  inner  cube  using  .25"  aluminum  plate  and  a  .25" 
space  between  cubes.  The  need  for  the  3"  inner  cube  was 
established  by  the  size  of  the  readily  available  .25"  square 
FSR,  Interlink  Electronics™  product  #30-301.  Each  sensor  had 
an  overall  length  of  2".  In  order  to  mount  these  on  the  inner 
cube  and  allow  enough  room  for  sensor  protective  padding  and 
sensor  wiring,  the  3"  dimension  was  chosen. 

The  four  FSRs  were  configured  as  shown  in  Figure 
23,  allowing  for  the  detection  and  resolution  of  the  four 
elements  of  the  force  and  moment  vector  according  to  Equation 


3 .  The  sensors  were  mounted  to  the  faces  of  the  inner 
plexiglass  cube  using  Scotch^  467MP  High  Performance 
adhesive.  A  10/32  threaded  set  bolt  was  inserted  above  each 
sensor  site  to  allow  for  the  application  of  ar  initial 
compressive  force  as  a  pre-load. 

To  protect  the  sensing  surface  of  each  FSR,  a 
.5''xl.75''  steel  shim  was  fastened  to  the  inside  face  of  the 
outer  cube.  The  sensor  side  of  this  shim  was  initially  lined 
with  .25"  neoprene  padding  material.  This  padding  was  later 
changed  to  a  .30"  polyurethane  bumpon”  to  better  protect  and 
transmit  normal  forces  to  the  FSR.  This  alteration  provided 
enough  padding  in  the  gap  between  cubes  to  exert  a  sufficient 
pre-compressive  load  to  eliminate  the  set  bolts  except  as  a 
method  for  manually  controlling  or  refining  the  set  point, 
b.  Electronic  Interface 

The  electronic  interface  shown  in  Figure  12  was 
used  to  convert  the  change  in  resistance  provided  by  each 
sensor  to  a  suitable  voltage  change.  The  se.nsors  were  wired 
such  that  FSR  #1  used  orange/yellow  colored  wires,  FSR  #2  used 
brown/ red  colored  wires,  FSR  #6  used  blue/green  colored  wires, 
and  FSR  #12  used  gray /violet  colored  wires.  These  were  multi¬ 
strand,  24  gauge  wires  separated  from  16  wire  ribbon  cable. 
Various  feedback  resistance  values  were  used  in  conjunction 
with  an  LM324  op-amp.  This  circuit  was  essentially  four  of 
the  circuits  used  in  the  pre-load  sensor  test  apparatus. 
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3.  Calibration  and  Taatiag 
a.  J^paratus 

A  general  sketch  of  the  calibration  and  test 
apparatus  is  shown  in  Figure  26.  The  device  was  constructed 
using  recdily  available  pine  l^xlO"  naterial  with  .125"  press 


Figure  2* 


Transducer  Test  and  Calibration  Apparatus 

board  used  to  reinforce  the  corners  as  shown.  The  uprights 
were  8.5"  high  and  the  device  was  24"  wide.  The  force-torque 
transducer  was  mounted  as  shown  enabling  individual  loading  in 
each  direction  of  interest.  A  system  consisting  of  two 
aluminum  shafts  mounted  through  ball  bearing  supports  was 
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added  to  reduce  friction  between  the  uprights  and  load 
supporting  strings  which  pass  over  the  ends. 
b.  Technique 

Incremental  masses  hung  from  strings  laid  over  the 
upright  ends  and  attached  to  the  joystick  handle  served  to 
load  the  system  in  the  x  and  y  directions.  The  joystick  cube 
had  be  rotated  90°  for  measuring  the  positive  and  negative  y 
forces  after  the  positive  and  negative  x  direction  forces  were 
applied  and  measured. 

To  measure  forces  in  the  positive  z  direction 
incremental  masses  were  placed  in  a  tray  which  was  mounted  to 
the  top  of  the  joystick  handle.  To  measure  forces  in  the 
negative  z  direction,  the  apparatus  was  designed  to  be  turned 
over  across  a  gap  between  two  adjacent  parallel  tables.  The 
loads  were  applied  by  incrementally  plt^cing  masses  in  a  tray 
hanging  from  the  joystick  handle. 

A  bar  was  bolted  to  the  top  of  the  joystick  handle 
parallel  to  the  upright  ends  in  order  to  measure  moments  about 
the  z  axis.  Masses  were  hung  in  opposite  directions  from  each 
end  of  this  bar  using  strings  passed  over  the  rollers  on  the 
uprights  providing  a  couple  with  the  bar  length  as  the  moment 
arm.  The  right  hand  rule  was  used  to  determine  the  directions 
in  which  to  load  the  bar  for  the  respective  positive  and 
negative  z  axis  moments. 
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c.  Data  Acquisition 

Masses  ranging  from  -2.762  kg  to  -4-2.762  kg  in 
increments  of  .25  kg  were  used  to  calibrate  the  transducer. 
The  moment  arm  provided  by  the  bar  bolted  to  the  top  of  the 
joystick  handle  was  .0984  meters.  The  masses  were  carefully 
used  to  apply  isolated  forces  in  x,  y,  and  z  directions  as 
well  as  isolated  moments  about  the  z  axis.  For  calibration, 
the  OMEGA  BENCH™  data  acquisition  system  and  APPLE  MACINTOSH™ 
microcomputer  were  configured  as  a  voltmeter  to  detect  the 
voltage  response  for  each  sensor  as  each  load  in  each 
direction  of  interest  was  applied.  The  changes  in  voltage  for 
each  change  in  load  were  plotted  as  16  separate  curves.  The 
slopes  of  these  curves  were  used  to  relate  the  force  and 
moment  column  vector  to  the  sensor  output  according  to 
Equation  24  as  the  16  elements  of  the  4x4  A  matrix. 

With  the  elements  of  the  A  matrix  determined, 
Equation  13  could  be  used  to  resolve  arbitrary  applied  reduced 
order  forces  and  moments.  The  OMEGA  BENCH™  data  acquisition 
icon  system  was  configured  as  shown  in  Figure  27  in  order  to 
perform  the  matrix  multiplication  in  real  time.  The  first 
column  of  icons  represents  the  incoming  analog  signals 
provided  by  the  sensor  outputs  in  response  to  an  applied  load 
on  the  joystick  handle.  The  second  column  of  icons  shows 
calculation  boxes  in  which  the  sensor  voltage  signals  were 
summed  with  the  initial  pre-load  voltages.  This  resulted  in 
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a  change  in  voltage  response  to  the  applied  load  which  could 
be  multiplied  by  the  elemencs  of  the  inverted  A  matrix.  The 
multiplication  was  done  in  the  third  and  fourth  columns  of 
icons  by  first  multiplying  each  change  in  sensor  voltage 
signal  by  the  inverted  A  matrix  element  in  the  corresponding 
matrix  column.  Then  adding  each  matrix  row  of  products 
provided  each  of  the  four  outputs.  Each  force  or  moment 
output  was  displayed  in  real  time  using  a  meter  and  a  chart 
icon. 

d.  Results 

Two  calibration  tests  were  conducted.  The  first 
used  a  feedback  resistor,  with  a  value  of  20  kohms  and 
the  second  used  10  kohms  for  feedback  resistance.  In  both 
tests,  it  was  possible  to  obtain  very  linear  force  and  moment 
versus  voltage  curves,  although  the  curves  appeared  to  be 
somewhat  discontinuous  at  the  zero  points.  This  was  due  to 
the  physical  variations  in  first  measuring  the  response  to 
applied  loads  in  the  positive  directions  and  then  unloading 
the  axis  under  consideration  and  applying  the  loads  in  the 
negative  direction  in  order  to  measure  the  response.  This 
discontinuity  was  accounted  for  by  determining  the  slopes 
using  four  different  techniques.  The  first  was  slopes 
determined  using  a  least-squares  regression  line  over  the 
entire  range.  The  second  and  third  were  using  the  slopes  from 
just  the  positive  and  negative  portions  of  the  curves. 
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Figur*  27 

Icon  ArrangoBont  for  Real  Tiaa  Matrix  Multiplication 


60 


respectively.  The  fourth  was  taking  the  average  of  the  two 
slopes  determined  for  each  curve  using  the  second  and  third 
technique  of  isolating  the  respective  positive  and  negative 
ranges . 

Once  calibrated  using  the  20  kohm  feedback 
resistance,  the  procedure  discussed  in  the  previous  section, 
and  the  slopes  resulting  from  the  least-squares  regression 
line,  the  system  was  tested  by  applying  a  knomi  force  of  10  N 
in  the  positive  x  direction.  The  first  time  this  load  was 
applied  the  response  was  good.  That  is  the  display  of  forces 
in  the  y  and  z  directions  remained  less  than  about  2  N  and  the 
moment  display  registered  very  little  change  while  the  x 
direction  force  display  was  about  8.5  N  indicating  successful 
decoupling  of  the  forces  and  moments.  However,  the  system  did 
not  prove  to  be  repeatable.  When  a  force  of  10  N  was  applied 
in  the  negative  x  direction  the  force  and  moment  displays  in 
the  y  and  z  directions  began  to  show  a  substantial  response 
indicating  poor  force  and  moment  resolution.  The  inverted  A 
matrix  elements  resulting  from  the  various  slope  determination 
techniques  were  all  tried  in  the  data  acquisition  system  in  an 
attempt  to  produce  a  more  repeatable  arrangement.  None  seemed 
to  provide  a  satisfactory  result.  The  matrix  elements  from 
the  technique  of  determining  the  slope  from  only  the  positive 
load  application  outputs  seemed  to  qualitatively  provide  the 
best  force  and  moment  resolution. 
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The  calibration  procedure  was  repeated  using  the  10 
kohm  feedback  resistance  with  no  pre-load  set  bolts  in  an 
attempt  to  refine  the  data  taking  technique  and  overcome  this 
non-repeatability.  The  results  were  virtually  the  same. 

B.  2  INCH  RBOUCBO  ORDER  VERSION 

1.  Objective 

In  parallel  to  building  a  4"  reduced  order  force- 
torque  transducer,  it  was  desired  to  design  a  smaller  reduced 
order  force-torque  transducer  to  meet  the  specific  payload 
requirements  of  the  PUMA  560  robotic  manipulator.  This 
required  a  device  with  a  2"  outside  dimension  weighing  less 
than  5  lbs.  In  order  to  do  this,  a  shorter  force  sensing 
resistor  was  needed. 

2 .  Shortened  rsR 
a.  Methods 

The  smallest  FSRs  available  measured  about  1.5"  in 
length  overall.  It  was  desired  to  use  the  same  type  of  FSR 
used  in  the  4"  transducer  design  in  the  smaller  transducer 
since  these  provided  a  good  response  for  a  given  force 
application.  Since  the  outside  dimensions  of  the  transducer 
had  to  be  approximately  2",  the  inside  cube  could  only  be 
about  1.25”  to  allow  a  .1875"  gap  using  .1875"  aluminum  plate. 
Hence,  the  FSR  resistor  used  could  not  be  any  longer  than  1" 
to  allow  enough  clearance  for  wire  attached  to  the  terminals. 
Also,  the  FSRs  could  be  no  wider  than  .5". 
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The  .25"  square  FSRs  used  in  the  4"  transducer  had 
an  overall  length  of  2".  In  order  to  reduce  this  length,  two 
methods  were  tried  which  did  not  use  soldering  due  to  the  low 
melting  temperatures  of  the  materials  used  to  manufacture  the 
FSRs.  Both  methods  attempted  to  leave  the  sensing  pad  of  the 
device  undisturbed  and  only  shorten  the  tail  portion  which  was 
used  to  connect  to  the  circuit  wiring.  The  first  method 
consisted  of  separating  the  outside  polymer  layers  on  the  tail 
portion  of  the  device.  With  the  silver  electrode  terminals 
exposed,  a  1"  section  was  removed.  The  remaining  .25"  portion 
crimped  with  the  wire  terminals  was  inverted  and  attached  to 
the  exposed  electrodes  near  the  sensing  portion  with  a 
electrically  conductive  epoxy.  After  a  several  attempts  an 
adequate  bond  between  the  shortened  pieces  was  not  achieved. 

The  second  technique  again  required  splitting  the 
polymer  covering  layers  on  the  tail  portion  of  the  FSR.  The 
device  was  trimmed  such  that  only  .125"  of  exposed  electrode 
remained  as  the  tail  piece.  A  small,  .0625"x.l25"  piece  of 
the  polymer  was  trimmed  away  from  each  side  of  each  exposed 
silver  electrode.  Small  metal  terminal  pin  connectors  with  24 
gauge  multi-strand  wires  pre-soldered  in  position  were  then 
spread  with  conductive  epoxy  and  crimped  onto  the  exposed 
electrodes.  This  provided  a  strong  bond  and  wire  attachment 
for  the  shortened  FSRs. 
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b.  Teat  and  Compariaon 

To  determine  the  effects  of  this  shortening 
process,  a  shortened  FSR  was  placed  in  the  pre-load  sensor 
test  apparatus  shown  in  Figure  11  using  the  electronic 
interface  shown  in  Figure  12.  The  test  conducted  was  similar 
to  the  previous  test  conducted  to  prove  the  pre-load  concept. 
A  plot  of  the  results  is  provided  in  Figure  28.  The  results 
showed  a  very  linear  relationship  between  the  applied  force 
and  the  voltage  response  when  compared  to  the  earlier  pre-load 
concept  test  results  shown  in  Figure  14. 
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Shortened  FSR  Performance  Test  Results 
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It  was  observed  that  for  the  shortened  test  the 
time  between  incremental  loadings  was  relatively  short  and 
consistent.  During  the  pre-load  concept  test,  no  specific 
care  was  taken  to  provide  a  consistent  time  interval  between 
load  applications.  This  was  thought  to  be  important  since  at 
this  point  it  was  noticed  that  the  output  slowly  crept  higher 
as  the  device  was  loaded  and  left  undisturbed. 

The  test  was  repeated  using  a  feedback  resistance 
of  10  kohms  for  three  different  initial  pre-load  set  points. 
A  plot  of  the  results  is  provided  in  Figure  29.  Run  #2  and 


run  #3  appear  to  provide  the  best  set  point  since  they  offer 
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the  largest  change  in  voltage  for  a  given  change  in  load  and 
remain  roughly  linear.  At  this  point,  the  shorter  FSR  appears 
to  provide  a  more  linear  output  than  the  unaltered  sensor. 
This  was  initially  thought  to  be  the  result  of  consciously 
applying  the  incremental  loads  and  recording  the  voltage 
response  data  in  a  quick  consistent  fashion  in  order  to 
compensate  for  the  drift  phenomenon  noticed  earlier. 

To  test  the  effects  of  this  data  taking  technique, 
an  unaltered  FSR  was  placed  in  the  apparatus  and  loaded  to 
attain  three  set  points  similar  to  the  previous  test  on  the 
shortened  FSR.  Figure  30  shows  the  results  of  this  test.  The 
response  curves  are  nearly  linear  indicating  that  the  speed 
and  consistency  of  the  data  taking  technique  was  important  due 
to  the  drift  tendency.  The  important  difference  was.  that  the 
pre-load  set  points  used  in  the  shortened  FSR  test  were  not 
repeatable  in  the  unaltered  FSR  test  using  a  feedback 
resistance  of  10  kohms.  The  higher  set  point  was  attained 
using  a  feedback  resistance  of  30  kohms  while  the  lower  two 
set  points  were  made  using  a  20  kohm  feedback  resistance.  This 
indicated  that  the  shortening  process  added  about  10  to  20 
kohms  of  resistance  in  series  with  the  shortened  FSR  since  the 
unaltered  sensor  required  additional  resistance  to  attain  the 
same  three  set  points  used  earlier  and  shown  in  Figure  29. 
The  resistance  added  to  the  altered  FSR  is  thought  to  be  the 
result  of  the  conductive  epoxy  used  in  the  fastening  process 
of  the  wire  terminals. 
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Figure  30 


unaltorod  F8R  Tost  at  Tliraa  Diffaraat  Pra>Loads 

c.  Drift  Analysis 
(1)  Description 

As  discussed  in  the  previous  section,  the  drift 
phenomenon  was  noticed  while  applying  loads  and  recording  the 
voltage  response  during  testing  of  the  shortened  FSR.  The 
presence  of  voltage  response  drift  was  confirmed  by  leaving 
the  apparatus  loaded  and  undisturbed  for  a  period  of  about  3 
hours.  In  this  time  the  voltage  output  ranged  from  5.988 
volts  to  6.290  volts.  However,  when  an  arbitrary  compressive 
force  was  suddenly  applied  and  released,  the  voltage  output 
dropped  back  to  approximately  5.9  volts.  Therefore,  this 
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indicates  that  the  effects  of  drift  do  affect  the  output  and 
the  response  curves  whose  slopes  are  used  to  determine  the 
elements  of  the  K  matrix. 

(2)  Test  Conducted 

To  more  fully  understand  the  effects  of  drift, 
an  arrangement  consisting  of  the  four  sensors  of  the  4" 
prototype  transducer  with  a  feedback  resistance  of  10  kohms, 
an  unaltered  FSR  in  the  pre-load  test  apparatus  with  a 
feedback  resistance  of  30  kohms,  an  unaltered  FSR  with  no  load 
with  a  feedback  resistance  of  l  Mohms,  a  shortened  FSR  loaded 
with  vice  grips  using  minimal  padding  with  a  feedback 
resistance  of  1  Mohms,  and  an  open  channel  with  5  volts 
directly  applied  using  no  FSR  with  a  feedback  resistance  of  l 
Mohms  was  developed.  Therefore,  this  arrangement  consisted  of 
eight  channels,  numbered  channel  1  to  channel  8,  representing 
the  eight  respective  sensor  setups  just  delineated.  The 
voltage  output  was  plotted  over  a  specific  period  of  time  at 
various  intervals. 

Two  tests  were  conducted  using  this  eight 
sensor  arrangement.  The  first  test  was  over  a  period  of 
almost  65  hours  with  17  output  recordings  at  various  time 
intervals  for  all  eight  channels  with  little  delay  between 
channels.  During  this  test  a  sudden  compressive  load  was 
applied  and  released  twice  and  a  facility  power  outage 
occurred.  The  second  test  covered  a  175  minute  period  with 
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the  output  data  taken  regularly  every  10  elnutes  for  each  of 
the  respective  eight  channels. 

Another  test  was  conducted  using  a  single 
unaltered  .25"  square  FSR  in  order  to  further  tiitudy  the 
effects  of  the  padding  used  between  the  load  application  and 
the  FSR  active  surface.  This  test  consisted  of  placing  the 
FSR  between  two  smooth  rigid  metal  surfaces.  One  was  aluminum 
plate  and  the  other  was  steel  plate  material.  The  data  was 
taken  first  over  a  period  of  276  minutes  in  13  intervals. 
These  were  every  3  minutes  for  the  first  27  minutes  and  then 
at  the  246,  261  and  276  minute  points.  A  second  set  of 
readings  was  taken  over  a  70  minute  period  in  23  Interval. 
These  intervals  were  every  minute  for  the  first  10  minutes  and 
increased  to  every  5  minutes  for  the  remainder  of  the  time. 

(3)  Results 

The  results  of  the  four  tests  discussed  in  the 
previous  section  are  provided  in  Appendix  D.  The  first  test 
was  conducted  over  the  longest  interval.  The  general 
conclusions  from  the  data  plots  is  that  the  drift  is  apparent 
in  all  channels  except  channels  #6  and  /8  which  were  the 
channels  with  the  unloaded  FSR  and  no  FSR,  respectively.  It 
was  also  apparent  that  as  the  load  condition  was  disturbed,  by 
applying  a  sudden  change  and  quickly  returning  the  original 
loading  condition,  the  voltage  value  would  drop  back  to  the 
initial  value  and  resume  drifting  upward.  This  load  change 
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was  experienced  by  the  apparatus  at  elapsed  tines  of  10.37 
hours  and  57.11  hours.  Also,  at  59.46  elapsed  hours  the 
facility  experienced  a  power  outage  which  caused  the  voltage 
value  to  drop  and  resume  the  upward  drift.  Throughout  all 
these  disturbances,  the  outputs  of  channels  #6  and  #8  remained 
nearly  constant. 

The  second  test  conducted  was  similar  to  the 
first  except  that  the  overall  time  period  was  shorter  and  the 
interval  between  voltage  readings  was  more  regular  at  10 
minutes  each  except  toward  the  end  of  the  test.  The  results 
were  as  expected,  demonstrating  the  drift  in  all  but  channels 
#6  and  #8.  The  only  disturbance  experienced  was  an 
intentional  system  power  outage  after  155  minutes  had  elapsed. 
It  was  at  this  point  that  the  data  acquisition  interval  was 
decreased  to  5  minutes.  Except  for  channels  #6  aiid  #8  the 
voltage  outputs  show  the  drift  upward  both  after  initial  load 
application  and  after  the  power  disturbance.  In  the  resulting 
plots  the  approach  of  the  effected  output  curves  can  be  seen 
to  approach  some  apparent  asymptotic  voltage  level. 

In  order  to  understand  the  effects  of  the 
padding  between  the  applied  load  and  the  FSR  sensing  surface, 
the  last  two  tests  were  conducted  as  discussed  in  the  previous 
section.  The  results  of  the  tests  clearly  indicate  a  voltage 
output  gradually  rising  toward  some  asymptotic  value.  No 
disturbances  were  caused  or  occurred  during  these  two  tests. 
The  conclusion  was  therefore  drawn  that  the  type  of  padding 
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between  the  load  and  the  FSR  sensing  surface  had  no  effect  on 
the  upward  drift  experienced  by  the  voltage  output. 

In  the  eight  channel  test,  Channel  #8  was  used 
to  show  that  this  drift  phenomenon  was  not  caused  by  the 
circuitry  used  since  no  drift  was  experienced  when  no  FSR  was 
used.  The  conclusion  drawn  from  the  results  of  the  Channel  #6 
test  is  that  an  unloaded  FSR  does  not  experience  a  drift  in 
response.  Hence,  it  was  shown  that  only  FSRs  with  applied 
loads  experience  a  circuit  output  drift  upward  toward  some 
final  steady  state  voltage  value.  This  drift  can  be  stopped 
and  started  again  from  an  initial  voltage  output  value  by 
suddenly  disturbing  the  loading  or  interrupting  the  system 
power . 

3.  Design  and  Construction 
a.  Mechanical  Device 

The  design  of  the  2”  force-torque  transducer  was 
very  similar  to  that  of  the  4”  force-torque  transducer.  The 
detailed  mechanical  drawings  are  provided  in  Appendix  E.  The 
primary  difference  other  than  the  size  was  the  round  top  and 
bottom  to  allow  enclosure  of  the  device  in  a  cylindrical  tube. 
Also,  the  inner  cube  was  made  from  aluminum  with  an  aluminum 
joystick  handle  threaded  and  inserted  into  a  threaded  hole  in 
the  top.  The  top  of  the  handle  was  threaded  to  allow  for 
fastening  the  bar  used  to  measure  the  moments  about  the  z 
axis.  The  size  of  the  inner  cube  allowed  installation  of  the 


shortened  FSRs  exactly  on  the  edges  in  accordance  with  the 
reduced  order  configuration  shown  in  Figure  24.  This  layout 
was  chosen  so  that  pre-load  set  bolts  could  be  applied  to  the 
sensor  in  the  z  direction  at  the  top  of  the  cube,  which  is  the 
negative  z  face,  in  order  to  allow  mounting  the  transducer  on 
the  end  effector  of  the  PUNA  560  robotic  manipulator. 
Polyurethane  bumpons*  which  were  .44'*x.20*'  were  used  between 
cubes  in  order  to  apply  normal  forces  to  the  sensors.  Since 
these  were  slightly  larger  than  the  .1875"  gap,  they  provided 
enough  pre-compress ive  load  after  completely  assembling  the 
device  to  eliminate  the  need  for  the  set  bolts.  The  set  bolts 
were  therefore  only  necessary  as  a  means  of  fine  tuning  the 
initial  set  points  by  adjusting  the  individual  pre-loads  on 
each  sensor. 

b.  SlBctronic  Inter facB 

The  current-to-voltage  circuit  shown  in  Figure  12 
was  again  used  for  this  smaller  design.  Feedback  resistance 
values  of  10  and  20  kohms  were  used.  The  circuit  design  was 
essentially  identical  to  that  used  for  the-  4"  prototype 
transducer  discussed  earlier  in  this  chapter.  The  only 
difference  was  the  wiring  color  scheme  used.  FSR  #3  used  the 
gray/ violet  colored  wires,  FSR  #4  used  the  white/black  colored 
wires,  FSR  #7  used  the  brown/red  colored  wires,  and  FSR  #10 
used  the  orange/yellow  colored  wires.  The  wires  were  the  same 
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type  of  24  gauge  multi-strand  insulated  wires  cut  from  16  wire 
ribbon  cable. 

4.  Calibration  and  Testing 

The  technique  of  calibrating  the  2**  transducer  was 
virtually  the  same  as  that  of  the  4"  transducer.  The 
apparatus  shown  in  Figure  26  and  the  OMEGA  BENCH™  icon 
arrangement  shown  in  Figure  27  were  designed  to  be  used  with 
both  the  2”  and  4"  devices  in  a  manner  similar  ro  that 
described  for  the  4"  version  earlier.  Four  calibration  and 
testing  attempts  were  made. 

The  first  two  calibration  tests  were  conducted  using 
a  feedback  resistance  value  of  10  kohms.  The  slopes  of  the 
resulting  16  least-squared  regression  lines  were  calculated 
for  each  of  the  two  test  runs.  The  respective  sets  of  16 
slopes  were  used  as  the  elements  of  two  A  matrices.  These 
matrices  were  inverted  and  used  in  the  calculation  icon  boxes 
to  perform  the  real  time  matrix  multiplication  in  accordance 
with  Equation  13.  The  resulting  output  was  similar  to  that  of 
the  4”  transducer  when  an  arbitrary  load  in  -the  x,  y  or  z 
direction  was  applied.  The  force  was  not  consistently 
resolvable.  As  the  joystick  was  manipulated,  force 
resolution,  if  attained,  was  not  repeatable.  As  arbitrary 
moments  about  the  z  axis  were  applied,  the  same  effects  were 
observed.  The  third  calibration  test  was  conducted  using  a 
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feedback  resistance  of  20  kohms.  The  results  were  also  not 


repeatable  with  no  clear  resolution  of  forces  or  moments. 

The  fourth  test  was  conducted  using  no  pre-load  set 
bolts  which  had  been  used  earlier  to  adjust  the  initial  pre¬ 
load  voltage  to  the  nominal  starting  value  used  during 
calibration.  This  test  showed  very  good  resolution  of 
arbitrarily  applied  forces  in  the  negative  direction,  but 
still  no  repeatable  results  were  attained. 
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VZZ.  DZSCnSSZOM  OF  RBBULT8 


A.  DRZFT  FHEMOMBMOM 

The  presence  of  drift  in  output  voltage  over  time  was 
demonstrated  in  the  previous  chapter.  This  phenomenon  was 
discovered  while  attempting  to  calibrate  the  2"  FSR  force-* 
torque  prototype  transducer.  The  drift  in  voltage  output  was 
noticed  while  loading  the  device.  The  output  crept  upward  as 
the  load  was  left  undisturbed.  The  output  would  drop  and 
resume  the  upward  drift  from  a  voltage  very  near  the  starting 
voltage  as  a  result  of  disturbing  the  load.  The  load 
disturbances  during  the  tests  consisted  of  suddenly  applying 
and  removing  a  larger  load.  Also,  an  interruption  in  power 
would  cause  the  drift  to  repeat.  In  all  cases  the  drift  in 
voltage  appeared  to  approach  some  asymptotic  value. 

The  effects  of  this  drift  were  not  accounted  for  in  the 
initial  test  of  the  pre-load  concept.  The  output  is  shown  in 
Figure  14  and  does  not  appear  as  linear  as  the  output  shown  in 
Figure  28.  The  output  shown  in  this  later  figure  was  produced 
by  attempting  to  compensate  for  the  drift  by  using  equal 
increments  of  mass  additions  over  roughly  equal  intervals  of 
time.  In  an  effort  to  remain  consistent,  the  emphasis  was  on 
preventing  any  disturbance  which  would  disrupt  the  drift  while 
taking  data.  The  artificiality  imposed  by  this  rigid  data 
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talcing  technique  was  duplicated  during  the  calibration 
procedures  conducted  on  the  prototype  transducers.  However, 
In  actual  use  of  the  transducers,  It  would  not  be  possible  to 
accurately  account  for  the  drift  phenomenon.  Sudden  changes 
In  loading  are  Inherent  In  the  operation  of  force-torque 
transducers . 

Therefore,  It  was  determined  to  disregard  the  effects  of 
the  drift  phenomenon  while  developing  physical  models  of  the 
force-torque  transducers  based  on  force  sensing  resistor 
technology.  Whan  a  prototype  Is  built  which  Is  capable  of 
repeatedly  resolving  the  desired  forces  and  moments,  the  error 
resulting  from  the  drift  could  be  determined.  Not  only  would 
drift  effects  have  been  difficult  to  quantify,  during  the 
design  of  the  prototypes  It  was  considered  superfluous  since 
the  drift  appeared  to  be  less  than  lOt  of  the  pre-load  voltage 
set  points.  The  accuracy  of  the  FSRs  used  In  this  application 
were  unknown  and  hence  the  overall  uncertainty  of  the 
transducer  output  was  not  known.  At  this  point  In  the  design, 
accuracy  was  not  the  primary  concern.  The  goal  was  to  show 
actual  decoupling  of  the  elements  of  the  force  and  moment 
column  vector,  F,  according  to  Equation  13. 

B.  PROTOTYPE  PERFORNAilCE 

1.  Single  FSR  Prototype 

The  single  FSR  prototype  displayed  In  Figure  11  using 
the  electronic  Interface  shown  In  Figure  12  was  used  as  a  test 
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platform  for  the  pre-load  concept  discussed  in  Chapter  III  and 
for  the  shortened  FSR  discussed  in  Chapter  VI.  As  mentioned 
previously,  during  the  testing  for  the  single  pre-loaded  FSR 
care  was  not  taken  to  consistently  load  the  device  without 
disturbing  the  drift  inherent  in  the  voltage  output.  However, 
the  device  was  used  to  successfully  show  the  ability  of  the 
pre-loaded  FSRs  to  be  used  for  a  voltage  signal  change,  for 
not  only  compressive  loads  but  also  tensile  loads.  This 
result  was  significant  because  it  allowed  the  required  number 
of  sensors  in  the  transducer  design  to  be  halved.  This  was 
because  one  pre-loaded  sensor  was  now  used  to  measure  forces 
in  both  the  positive  and  negative  directions  of  the  axis  on 
which  it  was  mounted. 

The  second  use  of  the  single  FSR  prototype  was  that  of 
proving  the  successful  FSR  shortening  process.  A  total  of  10 
.25”  FSRs,  Interlink  Electronics™  product  #  30-301,  were 
shortened  as  discussed  in  the  previous  chapter.  After  each 
FSR  was  altered,  it  was  mounted  in  the  single  FSR  pre-load 
device  and  tested  to  ensure  a  adequate  response. 

2 .  Transducer  Prototypes 

The  results  of  the  single  FSR  test  device  were  applied 
to  the  three  dimensional  theoretical  and  physical  force-torque 
transducer  models.  The  impact  of  the  successful 
demonstrations  of  the  pre-load  concept  and  the  FSR  shortening 
process  were  significant.  Also,  the  effects  of  the  interface 
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between  the  force  application  and  the  FSR  sensing  surface 
proved  to  be  important. 

First,  the  pre-load  concept  was  extended  to  the 
theoretical  model  of  the  three  dimensional  transducer  allowing 
the  associated  A  matrix,  originally  introduced  in  Equation  2, 
to  be  square.  This  was  a  direct  result  of  halving  the  number 
of  sensors  required  to  detect  forces  in  each  direction.  It 
was  shown  for  the  full  order  case,  which  sought  the  forces  in 
each  of  the  x,  y,  and  z  directions  as  well  as  the  moments 
about  these  three  axes  that  the  force  and  moment  column 
vector,  F,  would  have  six  elements.  The  elements  could  be 
resolved  using  a  transducer  configuration  of  only  six  sensors. 
A  similar  demonstration  was  made  for  the  reduced  order  case 
which  sought  the  forces  in  the  three  directions  but  only  the 
moment  about  a  single  axis.  This  simplified  the  problem,  as 
shown  in  Equation  13,  by  allowing  a  square  A  matrix  which,  if 
constructed  so  as  to  provide  a  full  rank,  could  be  inverted 
and  use  to  solve  for  F  directly. 

Second,  the  ability  to  shorten  the  individual  FSRs 
provided  the  means  of  building  a  force-torque  transducer  of 
smaller  dimensions  to  meet  specific  robotic  manipulator 
payload  requirements.  The  smaller  design  was  analogous  to  the 
first  larger  physical  three  dimensional  design  except  that  it 
used  the  altered  FSRs  whose  performance  had  been  proven  on  the 
single  FSR  test  device.  The  FSRs  altered  in  the  manner 
discussed  in  Chapter  VI  gave  the  ability  to  proceed  with  the 
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design  of  an  adequate  force-torque  transducer  to  be  used  with 
the  PUNA  560  robotic  manipulator  without  the  time  and  expense 
of  ordering  custom  FSRs  manufactured  to  meet  the  specific  size 
requirements . 

Finally,  in  the  actual  construction,  calibration  and 
testing  of  the  2"  and  4"  prototype  transducers,  output 
repeatability  was  a  problem.  Both  prototype  devices  provided 
an  initial  response  which  showed  definite  decoupling  of  the 
forces  and  moments  of  the  F  vector  according  to  the  reduced 
order  theory.  However,  a  consistently  repeatable  response  was 
not  obtained  regardless  of  many  attempts  to  accurately 
determine  the  elements  of  the  A  matrix  during  calibration 
procedures.  Ultimately,  this  lack  of  repeatability  was 
attributed  to  the  tendency  of  the  inner  cube  to  become  skewed 
in  an  orientation  with  respect  to  the  outer  cube  which  was 
somewhat  different  than  the  original  orientation.  The  padding 
used  to  transfer  the  load  to  the  sensor  was  a  set  of 
polyurethane  hemispherical  devices  which  theoretically 
provided  a  single  point  of  contact  on  the-  FSRs  sensing 
surface.  In  actuality,  each  pad  was  compressed  a  small  amount 
providing  a  relatively  small  flat  spot  that  remained  in 
contact  with  the  sensing  surface.  In  either  case,  the  load 
across  the  FSR  was  applied  in  a  relatively  small  local  pattern 
which  was  free  to  move  across  the  FSR  as  the  joystick  handle 
was  moved.  Moving  the  joystick  caused  the  inner  cube  to  move 
with  respect  to  the  outer  cube.  The  padding  was  firmly 
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attached  to  the  outer  cube  while  the  FSR  sensing  surface  was 
securely  attached  to  the  inner  cube.  The  result  was  a 
situation  in  which  the  application  point  of  the  load  was  free 
to  change;  therefore,  the  output  of  the  FSR  response  was 
reflecting  a  change  in  the  force  application  footprint.  This 
force  application  footprint  was  different  than  experienced  in 
calibration  of  the  device  requiring  a  different  A  matrix  to 
resolve  the  force  and  moment  vector,  F. 

This  effect  could  be  proven  by  utilizing  a  different 
polyurethane  pad  above  each  FSR.  These  pads  should  not  be 
hemispherical,  but  cylindrical  in  shape  so  that  the  force 
application  footprint  on  each  FSR  surface  remains  constant 
over  the  relatively  small  range  of  motion  between  the  inner 
and  outer  cubes  of  the  transducer.  After  eliminating  the 
force  application  footprint  as  a  variable,  the  device  should 
be  calibrated  as  discussed  earlier  and  retested.  Upon 
successful  and  repeatable  resolution  to  the  force  and  moment 
vector,  the  force-torque  transducer  accuracy  could  be 
quantified.  The  device  would  then  be  ready  to  -be  placed  in  a 
force  override  rate  control  system  as  discussed  in  the 
preliminary  work. 

C.  REDHHDANCY  ALGORITHM 

The  FORTRAN  code  developed  in  Chapter  IV  demonstrates  the 
ability  of  a  full  order  force-torque  transducer  to  be  built 
with  multiple  levels  of  redundancy.  The  theoretical 
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development  of  the  redundancy  algorithm  could  also  be  applied 
to  a  reduced  order  transducer  requiring  moments  only  about  a 
single  axis  as  modelled  in  Chapter  VI.  Once  proven,  the 
physical  device  could  be  built  with  multiple  levels  of 
redundancy  by  simply  adding  more  FSRs  and  a  computer  interface 
which  would  conduct  the  search  for  useable  combinations  as  the 
need  arose.  The  only  limit  to  the  levels  of  redundancy  would 
be  imposed  by  the  space  limitations  dictating  the  possible 
number  of  sensors  in  the  physical  transducer. 

The  full  order  design  used  to  develop  the  algorithm 
allows,  at  a  minimum,  four  layers  of  redundancy  simply  by 
using  all  the  available  corners  of  the  cube.  In  the  sample 
c  tput  in  Appendix  B,  it  took  li  sensor  failures  to  exhaust 
the  four  initial  subsystems  set  up  by  using  24  FSRs  and 
requiring  six  sensors  for  each  subsystem.  This  left  13 
sensors  in  an  arbitrary  configuration  from  which  the  search 
for  subsequent  sets  of  six  sensors  would  be  conducted  to  find 
useable  configurations.  In  the  example  output,  seven  such 
useable  systems  were  found  in  less  than  one  minute  after 
eliminating  a  sensor  each  time  another  set  was  discovered. 
This  example  provided  for  18  sensor  failures  out  of  24 
possible  sensors.  The  maximum  levels  of  redundancy  for  this 
configuration  is  18  since  at  least  six  sensors  are  required  to 
resolve  the  six  elements  of  the  force  and  moment  vector,  P. 
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VZII.  COMCLU8IOM8  AMD  RBCOIOIBMDATZOM8 


A.  COMCLU8IOM8 


•  Force  Sensing  Resistors  offer  a  practical,  robust  and 
inexpensive  alternative  to  strain  gauge  technology  for  use 
in  force-torque  transducers. 

•  Two,  three  dimensional  reduced  order  force-torque 
transducers  were  successfully  designed,  constructed  and 
used  to  qualitatively  demonstrate  the  ability  to  resolve 
forces  and  moments. 

•  Three  dimensional  force-torque  transducers  are  capable  of 
multiple  levels  of  redundancy  using  a  simple  algorithm  to 
search  for  acceptable  configurations  given  an  adequate 
remaining  layout  of  sensors  after  a  sensor  fails. 

•  The  effects  of  voltage  output  drift  over  time  will  effect 
the  overall  accuracy  of  resolving  forces  and  moments  using 
FSR  based  force-torque  transducers. 

•  Effects  of  the  interface  between  the  force  application 
surface  and  force  sensing  resistor  significantly  impact 
the  repeatability  of  the  force  and  moment  resolution  for 
the  force-torque  transducer  due  to  position  instability 
preventing  a  consistent  force  application  footprint  on 
individual  sensors. 


B.  RECOM1CEMDATIOM8 

•  Conduct  a  study  of  different  force  application  surfaces 
and  quantify  the  effects  of  their  footprints  on  FSR  based 
force-torque  transducer  output  repeatability. 

•  Use  adequate  force  application  surface  in  existing  force- 
torque  transducer  design  to  obtain  a  repeatable  resolution 
of  the  element  of  the  force  and  moments  column  vector  of 
interest. 

•  Develop  application  software  to  utilize  the  resolved 
forces  and  moments  from  the  FSR  based  force-torque 
transducer  in  the  force  override  rate  control  system. 
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APPBMDIX  A 
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TEST  #1  OF  SINGLE  PRE-LOADED  PSRi 


RESISTANCE:  4/7  (kohma)  i  =  I  25 


DATA: 

MASS  (Ibm)  VOLTAGE  (v) 


1.97 

5.0 

1.49 

S.65 

1.73 

S.15 

1.47 

5.71 

.99 

S.26 

-.72 

6.53 

.51 

6.62 

.50 

6.63 

.48 

6.64 

.24 

6.74 

.24 

48 

7.06 

.51 

6.99 

.72 

7.11 

1.01 

7.16 

1.23 

7.20 

1.48 

7.26 

1.73 

7.29 

2.20 

7.38 

2.46 

7.45 

3.18 

7.52 

3.45 

7.61 

3.69 

7.62 

3.93 

7.65 

4.17 

7.69 

FORCE  (Ib£)i 


F  = 


32.n 


Response  Graph: 


(volts)  V. 


1 


3 


TEST  #2-#6  OF  SINGLE  PRE-LOAOEO  PSR  USING  VARIOUS  FEEDBACK 
RESISTANCES: 

RESISTANCE:  varioua  (kohnts)  i  =  1  14 


DATA: 

MASS  (9)  VOLTAGE  (v) 


-1149.8 
r  1043.1 
936.2 

724.1 

494.2 
-49.5 

.01 

.01 

49.5 

494.2 

724.1 

936.2 
1043.1 
1149.8 


VI. 

1 

V2  =V3  =V4  = 
1  1  1 

V5  = 
1 

20 

mpaaa 

23 

20 

20 

23  2Dl  23 

23 

apii 

^02323 

gni 

5.56 

2E2i 

3R3 

5.61 

2020123 

20 

^2320 

5.81 

^21123 

33 

5.84 

29Z3iI&i 

33 

20 

202^23 

33 

23 

23[£Dpil 

33 

23 

20(^20 

33 

2apiza 

33 

m 

SSBIhI 

Nh«r«  th«  corrasponding 
rasiatancas  ara  (kohma); 

R1  =2.5 
R2  =5.0 
R3  =7,5 
R4  =10.0 
R5  =12.5 


Forca  ia: 


F.  =m.  g  l0  -  (N) 


Raaponaa  Grapha: 


2.5  kohms 


“  5.0  kohms 


R«spons*  Graphs  cont 
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PROGRAM  TMESIS 


THESIS 

NAVAL  ENGINEERING 
FSR  FORCE/TORQUE  TRANSDUCER 
REDUNDANCY  ANALYSIS 

—  CharlM  A.  Gunzel-III,  LT.  USN 


PARAMETER  ( N-6 , NRA-N , NCA-N , LDA-N , LDGINV-N , NR-1 ) 

REAL  A(24.6),A1(6,6).A2(6,6) , A3 (6, 6) , A4 (6, 6) , 

+  A5(«,6),PAC(6,6) ,T0L,AMACH,CIMV1,CINV2,GINV3, 
GINV4,GZNV5,RNO(1000)  ,RUNF 

INTEGER  COL, ROH«  RND, I , J,K, RANKl , RANK2 , RANK3 , RANK4 , 

*  RANKS. FSR(24) ,GRP1.GRP2.GRP3,GRP4,NI.N2, 
XPOS.XNEG.YPOS.YNEG.ZPOS.ZNEG 
CHARACTER  ANS*1 

OPEN  (20, FILE- ' THESIS . DAT^ , STATUS- ' NEW ' ) 

LOAD  THE  MATRIX 

Usa  transducar  configuration  ahown  in  Figure  21. 


c  * 
c  * 
c  * 
c  * 
c  * 
c  * 
c  * 
c  * 
c  * 


DATA  (A(1,I),I-1,6)  /7. 5, 0,0,0, -1,1/ 

DATA  (A(2,I),I-1,6)  /-l.S, 0, 0, 0, 1, -1/ 

DATA  (A(3,I),Z-1,6)  /0,7 . 5, 0, 1, 0, 1/ 

DATA  (A(4,I),1-1,6)  /O, -1.5, 0,-1, 0,-1/ 

DATA  (A(5,I) ,1-1,6)  /4 . 5, 4 . 5, 3 , 1, -I, 0/ 

DATA  (A(17,I),I-1,6)  /4. 5, 4. 5, 3,-1, 1,0/ 

DATA  (A(7,I),1-1,6)  /-7 . 5, 0, 0, 0, 1, -1/ 

DATA  (A(8,I),I-1,6)  /I. 5, 0, 0, 0, -1, 1/ 

DATA  (A(9,I),I-1,6)  /O, -7 . 5, 0, -1, 0, -1/ 

DATA  (A(10,I),I-1,6)  /O, 1.5, 0,1, 0,1/ 

DATA  (A(11,I) ,1-1,6)  /-4. 5, -4. 5, -3, -1,1,0/ 
DATA  (A(23,I),I-1,6)  /4. 5. 4. 5, -3, 1,-1, 0/ 

DATA  (A(13,I) ,1-1,6)  /7 . 5,0, 0, 0, -1, -1/ 

DATA  (A(14,I) ,1-1,6)  /-I. 5, 0, 0, 0, 1, 1/ 

DATA  (A(15,I) ,1-1,6)  /O, 7 . 5, 0, 1, 0, -1/ 

DATA  (A(16,I),I-1,6)  /O, -1. 5, 0, -1, 0, 1/ 

DATA  (A(6,I),I-1,6)  /-4 . 5, 4 . 5, 3 , 1 , I, 0/ 

DATA  (A(18,I),I-1,6)  /4 . 5, -4 . 5, 3, -I, -1, 0/ 

DATA  (A(19,I),I-1,6)  /-7 . 5, 0, 0, 0, 1, 1/ 

DATA  (A(20,I) ,1-1,6)  /I. 5, 0,0, 0,-1,-!/ 

DATA  (A(21,I) ,1-1,6)  /O, -7. 5, 0, -1, 0, 1/ 

DATA  (A(22,I) ,1-1,6)  /O, 1 . 5, 0, 1, 0, -1/ 

DATA  (A(12,I),I-1,6)  /4 . 5, -4 . 5, -3 , -1, -1, 0/ 
DATA  (A(24,I) ,1-1,6)  /-4 . 5, 4 . 5, -3 , 1, 1, 0/ 

C  PRINT  THE  MATRIX 

15  WRITE  (*,*)'Tha  full  aatrix  ia  aa  follova: ' 
WRITE  (20,*)'Tha  full  aatrix  ia  as  follows:' 
WRITE  (*,*) 

WRITE  (20,*) 

WRITE  (*,400) ((A(ROW,COL) , COL-1, 6) , ROW-1, 24) 
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WRITE  (20,400) ((A(ROW,COL) , COL-1, 6) , ROW-1, 24) 

C  VERIFY  THE  MATRIX 

WRITE  (*,*) 

WRITE  (20,*) 

WRITE  (•,*)'Do  you  wish  to  chanqa  any  olanantB?' 

WRITE  (20,*)* Do  you  wiah  to  changa  any  •laBants?' 

READ  (*,300)  ANS 
WRITE  (20,350)  ANS 
IF  (ANS.EQ.'Y*)  THEN 
GO  TO  25 
ELSE 
GO  TO  35 

25  WRITE  (*,*)*Ent«r  the  row  and  column  number  of  the  element:* 
WRITE  (20,*)* Enter  the  row  and  column  number  of  the  element:' 
READ  (*,*)  ROW, COL 
WRITE  (20,*)  ROW, COL 
WRITE  (*,*) 

WRITE  (20,*) 

WRITE  (*,100) 

■f  'Enter  the  corrected  value  of  element  (*  ,ROW,  * ,  *  ,COL,  * ) : 

WRITE  (20,100) 

*  'Enter  the  corrected  value  of  element  ( * ,ROW, * , * ,C0L, * ) : 

READ  (*,*)  A(ROW,COL) 

WRITE  (20,*)  A(ROW,COL) 

WRITE  (*,*)*Do  you  wieh  to  change  another  element?* 

WRITE  (20,*) 'Do  you  wleh  to  change  another  element?* 

READ  (*,300)  ANS 
WRITE  (20,350)  ANS 
IP  (ANS.EQ.'Y*)  THEN 
GO  TO  25 
END  IP 
END  IP 

WRITE  (*,*)*Do  you  wish  to  view  the  matrix  again?* 

WRITE  (20,*)*Do  you  wish  to  view  the  matrix  again?* 

READ  (*,300)  ANS 
WRITE  (20,350)  ANS 
IP  (ANS.EQ.'YM  then 
GO  TO  15 
END  IP 

35  CONTINUE 

C  BREAK  PULL  MATRIX  INTO  4  SUB  MATRICES 

DO  37  ROW-1,6 
DO  36  COL-1,6 
A1 (ROW, COL) -A (ROW, COL) 

A2(ROW,COL)-A((ROW>6) ,COL) 

A3  (ROW,COL)-A(  (R0W-fl2)  ,COL) 

A4  (ROW,  COL)  -A (  (ROW-M8)  ,  COL) 

36  CONTINUE 

37  CONTINUE 

C  DISPLAY  THE  4  SUB-SYSTEMS 

WRITE  (*,*)*The  A-POSITIVE  face  system  is:* 

WRITE  (20,*) 'The  A-POSITIVE  face  system  is:* 

WRITE  (*,400) ((Al(ROW,COL) , COL-1, 6) , ROW-1, 6) 

WRITE  (20,400) ((Al(ROW, COL), COL-1, 6) , ROW-1, 6) 
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WRITE  (*,*) 

WRITE  (20,*) 

WRITE  (•.•)'Th«  A-NEGATIVE  tacm  aystaa  is:' 

WRITE  (20,*)'Ths  A-NEGATIVE  facs  systSB  is:' 

WRITE  (*,400) ((A2(ROW, COL), COL-1, 6) , ROW-1, 6) 

WRITE  (20,400) ( (A2 (ROW, COL) , COL-1, 6) , ROW-1, 6) 

WRITE  (*,*) 

WRITE  (20,*) 

WRITE  (*,*)'Ths  B-POSITIVE  facs  systss  is:' 

WRITE  (20,*)'Ths  B-POSITIVE  facs  systsB  is:' 

WRITE  (*,400) ((A3 (ROW, COL), COL-1, 6) ,  ROW-1, 6) 

WRITE  (20,400) ((A3 (ROW,COL), COL-1, 6), ROW-1, 6) 

WRITE  (*,*) 

WRITE  (20,*) 

WRITE  (*,*)'Ths  B-NEGATIVE  facs  systSB  is:' 

WRITE  (20,*)'Ths  B-NEGATIVE  facs  systSB  is:' 

WRITE  (*,400) ((A4(ROW, COL), COL-1, 6) , ROW-1, 6) 

WRITE  (20,400) ( (A4 (ROW, COL) , COL-1, 6) , ROW-1, 6) 

WRITE  (*,*) 

WRITE  (20,*) 

C  CHECK  THE  RANK  OP  EACH  OF  THESE  SUBSYSTEMS 
TOL-100*AMACH(4) 

CALL  LSGRR  (NRA,NCA,Al,LDA,TOL,RANKl,GINVAl,LDGINV) 
WRITE  (*,500)'Ths  rank  of  ths  A-POSITIVE  facs 
systSB  ia:',RANKl 

WRITE  (20,500) 'Ths  rank  of  ths  A-POSITIVE  facs 
systsB  is:',RANKl 

CALL  LSGRR  (NRA,NCA,A2,LDA,TOL,RANK2,GINVA2,LDGINV) 
WRITE  (*,500) 'Ths  rank  of  ths  A-NECATIVE  facs 
t  systSB  la:«,RAMK2 

WRITE  (20,500) 'Ths  rank  of  ths  A-NEGATIVE  facs 

*  systsB  is:',RANK2 

CALL  LSGRR  (NRA,NCA,A3,LDA,TOL,RANX3,GINVA3,LDGINV) 
WRITE  (*,500) 'Ths  rank  of  ths  B-POSITIVE  facs 
•t-  systsB  is:',RANK3 

WRITE  (20,500) 'Ths  rank  of  ths  B-POSITIVE  facs 

*  systSB  is:', RANKS 

CALL  LSGRR  (NRA,NCA,A4,LOA,TOL,RANK4,GINVA4,LOGINV) 
WRITE  (*,500) 'Ths  rank  of  ths  B-NEGATlVE  facs 

*  systSB  la:',RANK4 

WRITE  (20,500) 'Ths  rank  of  ths  B-NEGATIVE  facs 
•f  systsB  is:',RANK4 

IP  ((RANK1.EQ.6) .AND. (RANX2.EQ.6) .AND. 

+  (RANK3.EQ.6) .AND. (RANK4.EQ.6))THEN 
WRITE  (*,*) 

WRITE  (20,*) 

WRITE  (*,*)'Each  of  thess  subsystsBS  is  invsrtibls. ' 
WRITE  (20,*) 'Each  of  thsss  subsystsBs  is  invsrtibls.' 

END  IP 
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C  INITIALXZI  24  FSR  0*S  TO  ZERO 

DO  39  1-1,24 
PSR(X)-0 
39  CONTINUE 

C  GENERATE  A  RANDOM  FSR#  FROM  1  TO  24 

GRPl-0 

GRP2-0 

GRP3-0 

GRP4-0 

XPOS-4 

XNEG-4 

ypos-4 

YNEG-4 

ZPOS-4 

ZNEG-4 

J-0 

45 

F8R(J)-24*RNUMF() 

C  ENSURE  FSR#  0  Z8  NOT  A  CHOICE 

IF  (F8R(J) .BQ.O)  THEN 
J-J-1 
GO  TO  45 
END  IF 

C  ENSURE  THIS  IS  NOT  A  REPEAT  VALUE 

DO  46  I-l, (J-1) 

IF  (FSR(J) .EQ.FSR(I))  THEN 
J-J-1 
GO  TO  45 
END  IF 

46  CONTINUE 

C  REPORT  THE  FSR  CHOSEN  TO  FAIL 

WRITE  (*,*) 

WRITE  (20,*) 

WRITE  (*,500)  'Tha  sensor  I  will  sieulate  failed  is  FSR  #',FSR(J 
WRITE  (20,500) 'The  sensor  I  will  sieulate  failed  is  FSR  #',FSR(J 

C  TERMINATE  IF  THERE  ARE  NOT  ENOUGH  FSR'S  REMAINING 

IF  (J.EQ.19)  THEN 

WRITE (*,*)  'TERMINATION— >  Sorry,  there  are  not 
-f  enough  remaining  FSRs.' 

NRITE(20,*) 'TERMINATION— >  Sorry,  there  are  not 
enough  remaining  FSRs.' 

GO  TO  1000 
END  IF 

C  CHECK  WHICH  SUBSYSTEM  THIS  FAILED  FSR  IS  IN 

IF  ((FSR(J).GE.1).AND.(FSR(J).LE.6))  THEN 
GRPl-1 

WRITE(*,600)  'This  sensor  is  in  the  A-POSITIVE  face  system.' 
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WRITE(20, 600) 'This  ssnsor  Is  in  the  A-POSITIVE  face  systee. ' 
ELSE  IP  ((PSR(J).GE.7).ANO.(PSR(J).LE.12))  THEN 
GRP2-1 

WR1TE(*,600)  'Thie  sensor  is  in  the  A-NEGATIVE  face  systee.' 

WRITE(20, 600) 'This  seneor  xs  in  the  A*NEGATIVE  face  systee.' 

ELSE  IF  ((PSR(J) .GE.i3) .AND. (FSR(J) .LE.18))  THEN 
GRP3-1 

NRITE(*,600)  'Thie  eensor  is  in  the  B-POSITIVE  face  systee.' 

WRITE(20, 600) 'This  sensor  is  in  the  B-POSITIVE  face  systee.' 

ELSE 
GRP4-1 

WRITE(*,600)  'This  sensor  is  in  the  B'NEGATIVE  face  systee.' 

WRITE(20, 600) 'Thie  sensor  is  in  the  B'NEGATXVE  face  systee.' 

END  IF 

C  SET  THIS  FSR  SIGNAL  TO  ZERO 

DO  47  COL-1,6 
A(PSR(J) ,COL)-0.0 
47  CONTINUE 

C  CHECK  IF  THERE  ARE  ANY  COMPLETE  INTACT  SUBSETS  REMAINING 

WRITE  (*,*) 

WRITE  (20,*) 

IF  (GRPl.EQ.O)  THEN 

WRITE(*,*)  'The  A'POSITIVE  face  systee  is  still  active.' 

WRITE ( 20, *) 'The  A'POSITIVE  face  systee  is  still  active.' 

END  IF 

IF  (GRP2.EQ.0)  THEN 

WRITE(*,*)  'The  A-NEGATIVE  face  eystee  is  still  active.' 
WRITE(20,*) 'The  A-NEGATIVE  face  systee  is  still  active.' 

END  IF 

IF  (GRP3.EQ.0)  THEN 

WRITE(*,*)  'The  B-POSITIVE  face  systee  is  still  active.' 
HR1TE(30, *) 'The  B-POSITIVE  face  systee  is  still  active.' 

END  IF 

IF  (GRP4.EQ.0)  THEN 

HRITE(*,*)  'The  B-NEGATIVE  face  systee  is  still  active.' 
WRITE(20,*) 'The  B-NEGATIVE  face  systee  is  still  active.' 

END  IP 

IF( (GRPl.EQ.O) .OR. (GRP2.EQ.0) .OR. (GRP3.EQ.0) .OR. (GRP4 . EQ. 0) ) THEN 
GO  TO  55 
ELSE 

NR1TE(*,*)  'None  of  the  4  initial  subsystees  reeain  intact.' 
WRITE(20,*) 'None  of  the  4  initial  subsystees  reeain  intact.' 

GO  TO  56 
END  IF 

C  CHECK  IF  IT  IS  DESIRED  TO  TAKE  OUT  ANOTHER  FSR 

55  WRITE  (*,*) 

WRITE  (20,*) 

WRITE  (*,*)  'Press  1  if  you  wish  to  el ieinate  another  FSR?' 

WRITE  (20,*) 'Press  1  if  you  wish  to  elieinate  another  FSR?' 

READ  (*,*)  I 
WRITE  (20,*)  I 
IF  (I.EQ.l)  THEN 
WRITE  (*,*) 

WRITE  (20,*) 
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GO  TO  4S 
ELSE 

GO  TO  1000 
END  ir 

56  CONTINUE 

WNITE(*,*)'I  will  randoBly  pick  6  of  tho  r«Minin9  rsM  to 

*  gonorato  anothor  aubaystaa.  • 

IIRITE(20,*)  *Z  will  randoaly  pick  6  of  tha  raaaining  FSRa  to 
ganarata  anothar  aubayataa. ' 

WRITEC*,*) 

WRITE(20.a) 

C  CHECK  IP  THERE  REMAINS  AT  LEAST  ONE  PSR'S  IN  EACH  DIRECTION 
DO  57  I-1,J 

IP  ((PSR(J) .EQ.l) .OR. (PSR(J) .EQ.2) 

.OR.  (PSR(J}.BQ.13}.OR.  (PSR(J).EQ.14))  THEN 
XPOS-XPOS'l 
END  IP 

IP  ((PSR(J) .BQ.7) .OR. (PSR(J) .EQ.8) 

*  .OR.(PSR(J).EQ.19).OR.  (PSR(J).EQ.20))  THEN 

XNEC«XMBG-1 
END  IP 

IP  ((PSR(J) .EQ.3) .OR. (PSR(J) .BQ.4) 

*  .OR.  (PSR(J).BQ.15).OR.  (PSR(J).EQ.16))  THEN 

ypos-ypos-1 

END  IP 

IP  ( (PSR(J) .EQ.9) .OR. (PSR(J) .EQ.IO) 

*  .OR.  (PSR(J).EQ.21).OR.  (PSR(J).EQ.22))  THEN 

yNEG«yMBG-l 
END  IP 

IP  ((PSR(J) .EQ.S) .OR. (PSR(J) .EQ.6) 

*  .OR.  (PSR(J).BQ.17}.OR.(PSR(J).EQ.I8))  THEN 

ZPOS-ZPOS'l 
END  IP 

IP  (  (PSR(J)  .BQ.  11)  .OR.  (PSR(J)  .EQ.  12) 

*  .OR.  (PSR(J).BQ.23).OR.  (PSR(J).BQ.24))  THEN 

ZNEC-ZNEC*! 

END  IP 

57  CONTINUE 


IP  ((XPOS.EQ.O) .OR. (XNEG.EQ.O) 

.OR.  (yPOS.EQ.O)  .OR.  (yNEG.EQ.O) 

*  .OR.  (ZPOS.EQ.O)  .OR.  (ZNEG.EQ.O))  THEN 
WRITE  (*,*) 

WRITE  (20,*) 

WRITE  (*,*)  'TERMINATION— >  Thara  ara  not  anough  FSRa 
■f  raaaining  in  aach  dlraction. ' 

WRITE  (20,*) 'TERMINATION— >  Thara  ara  not  anough  FSRa 
■f  raaaining  in  aach  dlraction.' 

GO  TO  1000 
END  IP 

C  TERMINATE  IF  THERE  ARE  NOT  ENOUGH  PSR'S  REMAINING 

IP  (J.EQ.19)  THEM 
WRITE  (*,*) 

WRITE  (20,*) 

WRITE  (*,*)  'TERMINATION«»»>  Thara  ara  not  anough  FSRa 
raaaining  to  continua.' 
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WRITE  (20. *) *TERN1NAT10N>««>  Th«r«  ar«  not  onough  FSRs 

*  roMlnlng  to  continuo.' 

GO  :o  1000 
END  IP 

C  RANOONLy  PICK  «  OP  THE  REMAINING  PSR'S 

65  Nl-0 
N2-1 

66  CONTIN'JE 
XPOS-U 
XNEG-0 
VPOS-0 
YNBG-O 
ZPOS-0 
ZNIG-0 

DO  69  X-(J«1),24 

67  PSR(K)-24*RNUNP() 

C  ENSURE  TSSLi  0  IS  NOT  A  CHOICE 

IP  (PSR(X)  .BQ.O)  THEN 
GO  TO  67 
END  IP 

C  ENSURE  mis  IS  NOT  A  REPEAT  VALUE 
DO  66  I-1,J 

IP  ((PSR(I)).BQ.(PSR(X)))THEM 
GO  TO  67 
END  IP 

68  CONTINUE 

C  ENSURE  THERE  IS  ONLY  ONE  PSR  PER  PACE 

IP  ((P8R(X).EQ.l).OR. (PSR(X).BQ.2) 

*  .OR.  (PSR(X).BQ.13).OR.(PSR(X).EQ.14))  THEN 

XPOS-XPOSfl 
END  IP 

IP  ((PSR(X).EQ.7).OR.(PSR(X).EQ.8) 

•1^  .OR.(PSR(X).BQ.19).OR.  (PSR(X).EQ.20))  THEN 
XNEC«XNEG*»^1 
END  IP 

IP  ((PSR(X) .EQ.3) .OR. (PSR(X) .EQ.4) 

*  .OR. (PSR(X) .EQ.15} .OR. (PSR(X) .EQ.16))  THEN 

ypos-ypos-M 

END  IP 

IP  ((PSR(K) .BQ.9) .OR. (PSR(K) . EQ. 10) 

*  .OR. (PSR(X) .EQ.21) .OR. (PSR(K) . EQ. 22) )  THEN 

YNEG-yNEG4>l 
END  IP 

IP  ((PSR(X).EQ.5).OR. (PSR(X).EQ.6) 

*  .OR. (PSR(K) .EQ.17) .OR. (PSR(X) .EQ.18))  THEN 

ZPOS-ZPOS-fl 
END  IP 

IP  ((PSR(X) .EQ.ll) .OR. (PSR(K) .EQ.12) 
.OR.(PSR(X).EQ.23}.OR.  (PSR(X).EQ.24))  THEN 
ZNEG-ZNBG-fl 
END  IP 

IP  ((XPOSl.GT.I)  .OR.  (XNEGl.GT.l) 
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*  .OR. (YPOSl.GT.I) .OR. (YNEGl.GT.  1) 

*  .OR. (ZPOSl.GT.l) .OR.  (ZNEGl.GT.l))  THEN 

GO  TO  66 

ELSE  IF  (N1.EQ.6)  THEN 
GO  TO  70 
END  IP 

69  CONTINUE 

C  SET  UP  NEW  6X6  MATRIX 

70  CONTINUE 

DO  72  ROW«(J>l),K 
DO  71  COL-1,6 

A5((ROW-J) ,COL)-A(PSR(ROW) ,COL) 

71  CONTINUE 

72  CONTINUE 

C  CHECK  RANK  6  PRINT  MEM  MATRIX  IP  RANK-6 

CALL  LSGRR  (NRA,NCA, AS, LDA.TOL, RANKS, GINVAS.LDGINV) 

IP  (RANKS.lt. 6)  THEN 
N1-N141 

IP  (N1.EQ.N2*(10**4))THEN 
N2-N2-»^l 

tlRITE(*,800)  'Th«r«  hav«  b«*n',Nl,'  itaratlons  atteapting  to 
4  find  an  invertible  matrix  with  the  remaining  FSRs.' 

IIRITE(20,800)  'There  have  been'.Nl,'  iterations  atteapting  to 
4-  find  an  invertible  matrix  with  the  remaining  FSRs.' 

END  IP 
GO  TO  66 
END  IP 

WRITE  (*,*) 

WRITE  (20,*) 

WRITE  (*,800) 'After' ,Nl, '  iterations,  the  randomly 

*  generated  subsystem  is:' 

WRITE  (20,800) 'After', Nl,'  iterations,  the  randomly 

*  generated  subsystem  is: ' 

WRITE  (*,400) ((AS (ROW, COL) , COL-1, 6) , ROW-1, 6) 

WRITE  (20,400) ((AS (ROW, COL) , COL-1, 6) , ROW-1, 6) 

WRITE  (*,*) 

WRITE  (20,*) 

WRITE  (*,500) 'The  rank  of  this  new  subsystem  is:', RANKS 
WRITE  (20, 500) 'The  rank  of  this  new  subsystem  is:', RANKS 
WRITE  (*,*)  'This  subsystem  provides  an  invertible  matrix.' 
WRITE  (20,*) 'This  subsystem  provides  an  invertible  matrix.' 

C  CHECK  IP  IT  IS  DESIRED  TO  TAKE  OUT  ANOTHER  FSR 

WRITE  (*,*) 

WRITE  (20,*) 

WRITE  (*,*)  'Press  1  if  you  wish  to  eliminate  another  FSR?' 
WRITE  (20,*) 'Press  1  if  you  wish  to  eliminate  another  FSR?' 

READ  (*,*)  I 
WRITE  (20,*)  I 
IP  (I.EQ.l)  THEN 
J-J-fl 
XPOS-4 
XNEC-4 
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YP0S«4 

YNEC-4 

ZPOS-4 

ZNEG-4 

73  PSR(J)-34*RNUNF() 

C  ENSURE  PSR/  0  ZS  NOT  A  CHOICE 


IP  (PSR(J) .EQ.O)  THEN 
GO  TO  73 
END  IP 

C  ENSURE  THIS  IS  NOT  A  REPEAT  VALUE 

DO  74  I-l« (J-1) 

IP  ((PSR(I)) .EQ. (FSR(J)))THEN 
GO  TO  73 
END  IP 
74  CONTINUE 


60  TO  56 
ELSE 

60  TO  1000 
END  IP 


C  FORMAT  LIBRARY 


100 

200 

300 

350 

400 

500 

600 

700 

800 


FORMAT  (1X,A,I2,A,12,A) 
FORMAT  (1X«24(1X,F6.3)) 
FORMAT  (A) 

FORMAT  (10X,A) 

FORMAT  ((T18,6(F5.1,3X)}) 
FORMAT  (IX, A, 13) 

FORMAT  (IX, A, 12, A) 

FORMAT  (IX, 12) 

FORMAT  (IX, A, 17, A) 


1000  CONTINUE 
STOP 
END 
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Sample  Output... 
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Th«  full  Mtrlx  is  as  follows 


7.5  0.0  0.0  0.0  -1.0  1.0 

-1.5  0.0  0.0  0.0  1.0  -1.0 

0.0  7.5  0.0  1.0  0.0  1.0 

0.0  -1.5  0.0  -1.0  0.0  -1.0 

4.5  4.5  3.0  1.0  -1.0  0.0 

-4.5  4.5  3.0  1.0  1.0  0.0 

-7.5  0.0  0.0  0.0  1.0  -1.0 

1.5  0.0  0.0  0.0  -1.0  1.0 

0.0  -7.5  0.0  -1.0  0.0  -1.0 

0.0  1.5  0.0  1.0  0.0  1.0 

-4.5  -4.5  -3.0  -1.0  1.0  0.0 

4.5  -4.5  -3.0  -1.0  -1.0  0.0 

7.5  0.0  0.0  0.0  -1.0  -1.0 

-1.5  0.0  0.0  0.0  1.0  1.0 

0.0  7.5  0.0  1.0  0.0  -1.0 

0.0  -1.5  0.0  -1.0  0.0  1.0 

4.5  4.5  3.0  -1.0  1.0  0.0 

4.5  -4.5  3.0  -1.0  -1.0  0.0 

-7.5  0.0  0.0  0.0  1.0  1.0 

1.5  0.0  0.0  0.0  -1.0  -1.0 

0.0  -7.5  0.0  *-1.0  0.0  1.0 

0.0  1.5  0.0  1.0  0.0  -1.0 

4.5  4.5  -3.0  1.0  -1.0  0.0 

-4.5  4.5  -3.0  1.0  1.0  0.0 

Do  you  wish  to  chsngs  any  slsasnts? 

N 

Ths  A-POSITIVE  facs  systsa  is: 

7.5  0.0  0.0  0.0  -1.0  1.0 

-1.5  0.0  0.0  0.0  1.0  -1.0 

0.0  7.5  0.0  1.0  0.0  1.0 

0.0  -1.5  0.0  -1.0  0.0  -1.0 

4.5  4.5  3.0  1.0  -1.0  0.0 

-4.5  4.5  3.0  1.0  1.0  0.0 

The  A-NEGATIVE  face  systsa  is: 

-7.5  0.0  0.0  0.0  1.0  -1.0 

1.5  0.0  0.0  0.0  -1.0  1.0 

0.0  -7.5  0.0  -1.0  0.0  -1.0 

0.0  1.5  0.0  1.0  0.0  1.0 

-4.5  -4.5  -3.0  -1.0  1.0  0.0 

4.5  -4.5  -3.0  -1.0  -1.0  0.0 

The  B-POSITIVE  face  systsa  is: 

7.5  0.0  0.0  0.0  -1.0  -1.0 

-1.5  0.0  0.0  0.0  -  1.0  1.0 

0.0  7.5  0.0  1.0  0.0  -1.0 

0.0  -1.5  0.0  -1.0  0.0  1.0 

4.5  4.5  3.0  -1.0  1.0  0.0 

4.5  -4.5  3.0  -1.0  -1.0  0.0 

The  B-NE6ATIVE  face  systsa  is: 

-7.5  0.0  0.0  0.0  1.0  1.0 

1.5  0.0  0.0  0.0  -1.0  -1.0 

0.0  -7.5  0.0  -1.0  0.0  1.0 

0.0  1.5  0.0  1.0  0.0  -1.0 

4.5  4.5  -3.0  1.0  -1.0  0.0 

-4.5  4.5  -3.0  1.0  1.0  0.0 
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Th«  rank  of  tha  A-POSITIVE  faca  ayotaa  la:  6 

Tha  rank  of  tha  A-NEGATXVB  faca  aystaa  la:  6 

Tha  rank  of  tha  B-POSITZVE  faca  aystaa  is:  6 

Tha  rank  of  tha  B'NBGATZVE  faca  aystaa  is:  6 

Each  of  thasa  subsystaaa  is  Invartibla. 

Tha  sansor  I  will  siaulata  failad  is  FSR  /  4 

This  sansor  is  in  tha  A-POSITIVE  faca  aystaa. 

Tha  A'NECATZVE  faca  aystaa  is  still  activa. 
Tha  B-P08XTZVE  faca  ayataa  is  still  activa. 
Tha  B-NEGATXVE  fsca  aystaa  is  still  activa. 

Prass  1  if  you  wish  to  aliainata  anothar  PSR? 
1 


Tha  sansor  X  will  siaulata  failad  is  PSR  f  3 
This  sansor  is  in  tha  A-POSXTIVE  faca  aystaa. 

The  A-NE6ATXVE  face  aystaa  is  still  actiVa. 
Tha  B-POSXTXVE  faca  ayataa  is  still  activa. 
Tha  B'NEGATXVE  faca  aystaa  is  still  activa. 

Press  1  if  you  wish  to  aliainata  anothar  FSR? 
1 


Tha  sansor  X  will  siaulata  failad  is  FSR  F  18 
This  sansor  is  in  tha  B*POSXTIVE  faca  aystaa . 

Tha  A-NEGATXVE  faca  aystaa  ia  still  activa. 
Tha  B'NEGATIVE  face  systaa  is  still  activa. 

Prass  1  if  you  wish  to  aliainata  anothar  PSR? 
1 


The  sensor  I  will  siaulata  failed  is  FSR  #  22 
This  sensor  is  in  tha  B>NEGATIVE  faca  systaa. 

The  A-NEGATXVE  face  systaa  is  still  active. 

Press  1  if  you  wish  to  aliainata  anothar  FSR? 
1 


The  sensor  I  will  siaulata  failad  is  FSR  #  16 
This  sensor  is  in  tha  B-POSXTXVE  face  systaa. 

The  A-NEGATXVE  face  systaa  is  still  activa. 

Press  1  if  you  wish  to  aliainata  anothar  FSR? 
1 


The  sensor  X  will  siaulata  failad  is  FSR  #  21 
This  sensor  is  in  tha  B-NEGATIVE  faca  systaa. 
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Th«  A-NECATXVB  fac«  •ystM  is  still  sctivs. 


Prsss  1  if  you  wish  to  slisinsts  snothsr  PSR? 
1 


Ths  ssnsor  X  will  sisulsts  fsilsd  is  rSR  #  1 

This  ssnsor  is  in  ths  A-POSITIVE  fscs  systss. 

Ths  A'NBGATIVE  fscs  systsa  is  still  sctivs. 

Prsss  1  if  you  wish  to  sliainsts  snothsr  FSR? 
1 


Ths  ssnsor  X  will  sisulsts  fsilsd  is  PSR  #  20 
This  ssnsor  is  in  ths  B-NEGATIVE  fscs  systss. 

Ths  A-NEGATIVE  fscs  systsB  is  Still  sctivs. 

Prsss  1  if  you  wish  to  slisinsts  snothsr  PSR? 
1 


Ths  ssnsor  X  will  sisulsts  fsilsd  is  PSR  #  19 
This  ssnsor  is  in  ths  B-NEGATXVE  fscs  systss. 

Ths  A'NEGATXVE  fscs  systss  is  still  sctivs. 

Prsss  1  if  you  wish  to  slisinsts  snothsr  PSR? 
1 


Ths  ssnsor  I  will  sisulsts  fsilsd  is  PSR  #  2 

This  ssnsor  is  in  ths  A-POSXTXVE  fscs  systss. 

Ths  A-ME6ATXVE  fscs  systss  is  Still  sctivs. 

Prsss  1  if  you  wish  to  slisinsts  snothsr  PSR? 
1 


The  ssnsor  X  will  sisulsts  fsilsd  is  PSR  #  23 
This  isnsor  is  in  ths  B-NEGATXVE  fscs  systss. 

The  A-NEGATIVE  fscs  systss  is  still  sctivs. 

Prsss  1  If  you  wish  to  slisinsts  snothsr  FSR? 
1 


The  sensor  X  will  sisulste  fsiled  is  FSR  t  11 
This  sensor  is  in  the  A-NEGATXVE  fsce  systss. 

None  of  ths  4  initial  subsystess  resain  intact. 

X  will  randosly  pick  6  of  ths  rsaaining  PSRs  to  generate  another  subsystes. 


After  0  iterations,  ths  randosly  generated  subsystes  is: 
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4.5  7.5  1.0  1.0  0.0  -4.5 

1.5  0.0  -7.5  0.0  -1.0  0.0 

-4.5  -7.5  0.0  0.0  0.0  1.0 

7.5  4.5  4.5  3.0  1.0  -1.0 

0.0  7.5  0.0  0.0  0.0  -1.0 

1.5  1.5  0.0  0.0  0.0  -1.0 

Thtt  rank  of  thla  now  ouboyotoo  lo:  6 

This  auboyotoo  providoa  an  Invartibla  aatrix. 

Proaa  1  if  you  viah  to  olialnata  anothar  FSR? 

1 

I  will  randomly  pick  6  of  tha  raaaining  PSRa  to  ganarata  anothar  aubayataa 


Aftar  1  itarationa,  tha  randomly  ganaratad  aubayatam  ia: 

-7.5  0.0  -7.5  0.0  -1.0  0.0 

1.5  4.5  4.5  3.0  -1.0  1.0 

-7.5  -7.5  0.0  0.0  0.0  1.0 

0.0  -4.5  4.5  3.0  1.0  1.0 

-7.5  -4.5  4.5  3.0  1.0  1.0 

4.5  1.0  1.0  0.0  -4.5  1.0 

• 

Tha  rank  of  thia  naw  aubayataa  ia:  6 

Thia  aubayatam  providaa  an  invartibla  matrix. 

Praaa  1  if  you  wiah  to  aliminata  anothar  PSR? 

1 

I  will  randomly  pick  6  of  tha  ramaining  PSRa  to  ganarata  anothar  aubayatam 

Aftar  1  itarationa«  tha  randomly  ganaratad  aubayatam  ia: 

-4.5  4.5  -4.5  -3.0  -1.0  -1.0 

4.5  -1.5  0.0  0.0  0.0  1.0 

4.5  0.0  -7.5  0.0  -1.0  0.0 

-1.5  0.0  1.5  0.0  1.0  0.0 

-1.5  1.0  1.0  0.0  -4.5  1.0 

-4.5  4.5  3.0  1.0  1.0  0.0 

Tha  rank  of  thia  naw  aubayatam  ia:  6 

Thia  aubayatam  providaa  an  invartibla  matrix. 

Press  1  if  you  wish  to  eliminate  anothar  PSR? 

1 

1  will  randomly  pick  6  of  tha  raaaining  PSRa  to  generate  another  subsystem 

After  0  iterations,  tha  randomly  ganaratad  aubsystaa  is: 

4.5  0.0  1.5  0.0  1.0  0.0 

-4.5  4.5  4.5  3.0  1.0  -1.0 

4.5  4.5  4.5  3.0  1.0  -1.0 

-1.5  1.0  1.0  0.0  -4.5  1.0 

-4.5  4.5  3.0  1.0  1.0  0.0 

-1.5  0.0  0.0  0.0  1.0  1.0 

The  rank  of  this  naw  subsystem  is:  6 

This  subsystaa  provides  an  invartibla  matrix. 

Press  1  if  you  wish  to  aliminata  anothar  PSR? 

1 


I  vlll  randoaly  pick  6  of  tho  raaalnlng  FSRa  to  9anarata  another  subsyataa 


After  1  itoratlona,  tha  randomly  ganaratad  aubaystaa  la: 

4.5  0.0  >7.5  0.0  >1.0  0.0 

•1.5  4.5  4.5  3.0  1.0  >1.0 

4.5  1.0  1.0  0.0  -4.5  1.0 

4.5  4.5  3.0  -1.0  1.0  0.0 

0.0  -7.5  0.0  -1.0  0.0  -1.0 

1.5  0.0  0.0  0.0  -1.0  1.0 

Tha  rank  of  this  naw  aubaystaa  is:  6 

This  subsystsa  providsa  an  invartibla  matrix. 

Prsss  1  If  you  wish  to  sliainats  another  PSR7 
1 

I  will  randoaly  pick  6  of  tha  raaalning  PSRs  to  gsnarats  another  subsystaa 


After  1  iterations,  the  randomly  ganaratad  subsystaa  is: 

0.0  1.5  0.0  0.0  0.0  -1.0 

-4.5  1.0  1.0  0.0  -4.5  1.0 

-4.5  4.5  3.0  '  1.0  1.0  0.0 

1.5  0.0  0.0  0.0  -1.0  1.0 

4.5  4.5  3.0  -1.0  1.0  0.0 

-1.5  0.0  0.0  0.0  1.0  1.0 

The  rank  of  this  naw  subsystaa  is:  6 

This  subsystaa  provides  an  Invartibla  matrix. 

Press  1  if  you  wish  to  sliainats  another  FSR? 

1 

I  will  randoaly  pick  6  of  tha  remaining  FSRs  to  generate  another  subsystaa 


After  27  iterations,  tha  randoaly  ganaratad  subsystsa  is: 

-1.5  1.0  1.0  0.0  -4.5  1.0 

0.0  -7.5  0.0  -1.0  0.0  -1.0 

-1.5  0.0  0.0  0.0  1.0  1.0 

0.0  1.5  0.0  1.0  0.0  1.0 

4.5  4.5  3.0  -1.0  1.0  0.0 

1.5  0.0  0.0  0.0  -1.0  1.0 

The  rank  of  this  naw  subsystaa  is:  6 

This  subsystaa  provides  an  invartibla  matrix. 

Press  1  if  you  wish  to  sliainats  another  FSR7 
1 

I  will  randomly  pick  6  of  the  reaaining  FSRs  to  generate  another  subsystem 


TERMINATION«»»  There  are  not  enough  FSRs  remaining  to  continue. 
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APPENDIX  C 
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FORCE-TORQUE  TRANSDUCER 
(OUTSIDE  BC)X) 
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APPEMDZX  O 
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DRIFT  ANALYSIS  TEST  #1 


RESISTANCE:  various  (kohms)  i  -  I  17 


OATE/TIME  ELAPSED  TIME  (hrs) 


VOLTAGE  (V) 


dtg^  - 


vl  :  v2  :  v3  =  v4  :  v5  =  v6  :  v7  =  v8  “ 

I  I  I  I  1  1  II 


0 

1035 

10.37 

11.34 

12  44 

13.39 

14.76 

15.64 

1667 

23.10 

24.78 

57  08 

57.11 

58.23 

59  46 

62  34 

64.67 

6  149 

|6  090| 

6  161 

6  187 

6087 

004 

8704 

5059 

6  149 

6  176 

6220 

6216 

004 

8704 

5059 

5  749 

5975 

6089 

6089 

004 

8  701 

5059 

5963 

6.092 

6085 

6  170 

6  168 

004 

8709 

5060 

6009 

6  113 

6  115 

6.197 

6  193 

004 

8714 

5062 

6  041 

6  128 

6  140 

6.210 

6  211 

004 

8  719 

5062 

6060 

6  138 

6  164 

6238 

6230 

004 

8.722 

5  062 

6074 

6  142 

6  174 

6248 

6232 

004 

8724 

5062 

6.083 

6  147 

6  187 

6255 

6.240 

004 

8724 

5062 

6.062 

6  106 

6  182 

6.240 

6236 

004 

8709 

5064 

6060 

6  104 

6  178 

6234 

6.234 

004 

8704 

5064 

6011 

6050 

6  174 

6.197 

6253 

004 

8  691 

5060 

5.585 

5759 

5745 

5906 

6006 

004 

8.679 

5  060 

5  811 

5880 

5929 

6035 

6  109 

004 

8  681 

5060 

5852 

5  902 

5963 

6060 

6  138 

004 

8684 

5060 

5.844 

5887 

6048 

6  121 

004 

864(- 

5057 

5.923 

5950 

6031 

6.115 

6  176 

004 

868(< 

5060 

232247 

240908 

240909 

241007 

241113 

241210 

241332 

241425 

241527 

242153 

242334 

260752 

260754 

260901 

261015 

261308 

261527 

CHI— joystick  FSR#1,  lOkohms 

CH2— joystick  FSR#2, lOkohms 

CH3— joystick  FSR#12,  lOkohms 

CH4~  joystick  FSR#6,  lOkohms 

CH5 — uncut  FSR  in  beam  prototype, 30kohm8 

CH6 — uncut  FSR  w/  no  applied  load,lMohm8 

CH7 — 8hortened  FSR  loaded  with  vice  gripa,lMohm8 

CHS — no  FSR,  just  5v  directly  applied, IMohms 
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CHANNEL  1 


CHANNEL  2 


CHANNEL  7 


CHANNEL  8 


DATE /TIME 
dtg,  -- 


282010 

282020 

282030 

282040 

282050 

282100 

282110 

282120 

282130 

282140 

282150 

282200 

282210 

282220 

282230 

282240 

282245 

282250 

282255 

282300 

282305 

ELAPSED  TIME (min) 


t  ^ 

I 


DRIFT  ANALYSIS  TBST  #3t 


RESISTANCE:  1  (Mohms)  i  M  13 


DATE/TIME  ELAPSED  TIME(roin)  VOLTAGE  (v) 
dtg,  --  ‘i  =  v7.  = 


931611 

a 

931614 

a 

931617 

Bl 

931620 

931623 

m 

931626 

Q 

931629 

m 

931632 

m 

031635 

fSM 

931638 

m 

032017 

23 

032032 

032047 

276 

MMAfl 

im 

7.148 

7  168 

7.176 

7.187 

7.238 

7.263 

7.276 

7.498 

7.511 

7.513 
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